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Abstract 

In Differential Pulse Coded Modulation (DPCM) the 
prediction of the next sample value is formed from the past 
values. This prediction can be considered as a set of 
instruction for the quantizer to conduct its next sample 
value in a particular interval. By using the redundancy in 
the signal to form a prediction, the region of uncertainty is 
reduced and the quantization can be performed with a 
reduced amount of decisions (or bits) for a given 
quantization level or reduced amount of quantization levels 
for a given number of decisions (or bits). This paper 
basically focuses on the software modeling of an open-loop 
DPCM structure by using MATLAB. During Simulation, 
comparisons were made in terms of different Signal to 
Noise ratios by using different number of bits per sample 
of the uniform quantization process and the prediction 
orders and their optimal values were determined. By 
determining these optimal values of number of bits as well 
as the prediction orders, the implementation of open-loop 
DPCM in production scenarios such as lossless 
compression or embedded compression could be made 
much more efficiently. Standalone graphical user interface 
(GUI) was developed for this task. 1 

Keywords: dpcm, open-loop dpcm, quantization, 
differential pulse coded modulation, prediction error, pulse 
amplitude modulation, pulse coded modulation. 


Nomenclature 

DPCM Differential Pulse Coded Modulation 

m[k] K th sample of uniform quantization 

d[ k] Difference between two corresponding 

uniform quantization samples 
Av Quantization interval 

L Quantization level 

m[k\ Estimated value of K th sample 


1 This study has been implemented on MATLAB software 
at BG320/322 lab. La Trobe University. 


1. Introduction 

DPCM compression depends on prediction technique, well- 
conducted prediction techniques leads to good compression 
rates, in other cases DPCM could mean expansion 
comparing to regular PCM encoding. There are two 
structures for DPCM: open-loop and close-loop. This paper 
is mainly focused on the software implementation of the 
open-loop DPCM. 

DPCM is a computationally efficient means for waveform 
coding. Various common audio signals exhibit a long-term 
low-pass spectral characteristic which results in a 
significant correlation between successive samples of the 
input signal [1]. 

In order to address some of the applications of DPCM 
systems when it comes to audio compression, de-noising, 
video coding, prediction and communication, these papers 
could be considered [2, 3, 4, 5, 6]. 

Mansour et.al in their paper [7] created a simulation of the 
DPCM and ADM systems in Simulink. But there is still a 
lack of proper standalone GUI based software solutions to 
perform DPCM simulations. This research paper deals with 
this problem by creating a standalone GUI based open-loop 
DPCM simulation. Also this article deals with the 
determination of optimal values of bits per sample used in 
uniform quantization and the prediction orders by 
comparing SNRs. By doing so the applications mentioned 
in [2, 3, 4, 5, 6] could be performed more efficiently. 

Section 2 deals with the theoretical backgrounds applied to 
the simulation process. This includes defining the 
relationship between the successive samples and how 
transmitting the difference between them helps the overall 
system to become more efficient, rather than transmitting 
the signal itself. 

In section 3, the actual implementation of the simulation 
parameters were demonstrated. Various functions that were 
required to build up the robust GUI system for simulation 
were discussed in this section. Various combinations of 
experimental encoder and decoder values were also 
displayed in tabular format in this section. 

Section 4 deals with the simulation results and various plots 
were constructed on the basis of trade-offs between the bits 
per sample of the uniform quantization process and the 
prediction order in this section. Tables displaying the 
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various signal to noise ratios for trade-offs between the bits 
per sample of the uniform quantization process and the 
prediction order were displayed in this section such that 
their optimal values could be found out. 

Finally section 5 deals with the conclusion and possible 
scope of improvements for the project in future. 


2. Theoretical Fundamentals 

In analog messages it is possible to make a good guess about 
the sampled values from the knowledge of the past sample 
values. That is, the sample values are not independent and 
generally there is a great amount of redundancy in the 
Nyquist samples. Proper Exploration of this redundancy 
leads to encoding a signal with a lesser number of bits. If it 
is considered as a sample scheme where instead of 
transmitting the sample values, the system transmits the 
difference between the successive samples. Thus, if m[k] 
is the kth sample, instead of transmitting m[k\ , the system 
transmits the difference d[k\ = m[k\ — m[k — 1] [8]. 

At the receiver, knowing d[k] and the previous sample 
value m[k — 1], the system can reconstruct m[k]. So it can 
be see that from the knowledge of difference d[k\, the 
system can reconstruct m[k] iteratively at the receiver. Now 
the difference between the successive samples is generally 
much smaller than the sample values. Thus, the peak 
amplitude m p of the transmitted values is reduced 
considerably. Because the quantization interval Av =mp/L, 
for a given L (quantization level) this reduces the 
quantization interval Av, thus reducing the quantization 
noise, which can increase the SNR [9]. 

The proposed system can improve upon this existing 
scheme by estimating or predicting the value of the k th 
sample m[k] from a knowledge of the previous sample 
values. If this estimate is m[k\, then it actually transmit the 
difference (prediction error)d[/c] = m[k ] —m[k\. At the 
receiver also, we determine the estimate m[k ] from the 
previous sample values, and then generate m[k ] by adding 
the received d[k] to the estimate m[k\. So it is basically 
able to reconstruct the samples at the receiver iteratively. If 
everything has been done properly then the predicted or 
estimated value m[k ] will be close to m[k ] and their 
difference (prediction error) d [k] will be even smaller than 
the difference between the successive samples. 

3. Implementation 

In DPCM coding a difference is based on the fact that most 
source signals shows significant correlation between 
successive samples so encoding uses redundancy in sample 
values which implies lower bit rate. Realization of basic 
concept (described above) is based on technique in which 
we have to predict current sample value based upon 
previous samples (or sample) and the system has to encode 
the difference between actual value of sample and predicted 
value (difference between samples can be interpreted as 
prediction error). 


The open loop DPCM is comprised of mainly two main 
parts, the Encoder and the decoder. The encoding 
functionalities are implemented inside the encoder in the 
transmitting side while the decoding activities are 
performed by the decoder in the receiver side. Basically the 
encoder instead of encoding the input signal directly and 
sending it via the channel to the receiver, it encodes the 
difference between the main input signal sample and it’s 
linear prediction sample and sends it to the receiver for 
decoding. This is a much more efficient way of sending 
information from the transmitter and the receiver as it 
prevents a lot of redundant calculations. Due to this less 
number of bits will be needed to send the information via 
the channel so there will be less chance of bit error 
probability too. 

The functionalities inside the encoder can be explained with 
the help of following equations that occurs inside the 
encoder. 


M 


= 2^ a k x[n — k] 

( 1 ) 

k= 1 


e[n ] = x[n\ — x p [n] 

( 2 ) 

e q [n ] = Q(e[n]) 

( 3 ) 



Figure 1: Encoder (top) and Decoder (bottom) block diagrams 


Tablet : Encoder Experimental values 


n 

1 

2 

3 

4 

x[n] 

10 

20 

20 

40 

x P [n] 

0 

9 

18 

18 

e[n] 

10 

11 

2 

22 

e q [n] 

2 

2 

0 

5 


For the experimental implementation of the encoder, the 
system basically takes an input signal x[n\ = 
{10,20,20,40} The quantizer is designed as e q = 

floor(p-^). So it can be seen that after the input signal 
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traverses the entire block diagram of the encoder we are left 
with the values e q = (2, 2, 0,5}. The encoder basically 
transmit e q than the actual signal x[n\. 

The system can find the experimental values of the 
functioning of the decoder the same way as the encoder. 
The de- quantization block is designed as dQ = 4 * e q [n\. 
Following is the table showing the experimental values of 
the decoder block. 


Table 2: Decoder Experimental Values 


n 

1 

2 

3 

4 

e q [n] 

2 

2 

0 

5 

e[n] 

8 

8 

0 

20 


0 

7.2 

13.68 

12.312 

x[n\ 

8 

15.2 

13.68 

32.312 


The proposed system can also find the error due to 
excessive quantization from the above experimental values 
which is given by equation (4). For this the above case it is 
found to be E[n\ = (2, 5, 6, 7). 

E[n\ = x[n\ - x[n\ (4) 

Separate MATLAB functions entitled; 

• my Quantization, 

• myDeQuandzation, 

• myPrediction 

• myPredictionReconstrucdon 

were built to perform the functions of the block diagrams 
of the encoder and the decoder. A graphical user interface 
was also build that gave the user to choose between the 
three sound samples of handel, laughter and train in pop- 
up menu. 

At first the system loaded the sound file laughter with 
sampling frequency of 8192 Hz in the main mat file. The 
users were asked in the GUI to enter the Quantization level 
(K) as well as the Prediction Order (N). As advised in the 
handouts we performed the encoding and the decoding 
process by using the Prediction Order of 3 (N =3) and 
Quantization Level of 4 (K = 4). Following graphs were 
obtained: 



Figure 2: GUI for simulation of open-loop DPCM 


Original Signal 

1 1 1 1 1 — 



0 1000 2000 3000 4000 5000 6000 7000 

Time (ms) 

Figure 3: Original Signal (laughter) 

Quantized Signal 

1 1 — i 1 1 1 — i 1 



_1 .5 i 1 1 1 1 1 1 

0 1000 2000 3000 4000 5000 6000 7000 

Time(ms) 

Figure 4: Quantized Signal (laughter) 
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0.8 
0.6 
0.4 
0.2 
o 

- 02 . 

- 0.4 
- 0.6 
- 0.8 
-1 

0 1000 2000 3000 4000 5000 6000 7000 

Time(ms) 

Figure 5: Signal after De-quantization (laughter) 

Reconstructed signal by DPCM 


0.5 

S o 

- 0.5 

-1 

- 1 .5 

0 1000 2000 3000 4000 5000 6000 7000 

Time(ms) 

Figure 6: Signal after Linear prediction by order 3 (laughter) 

In a similar way the other two signals handel and train 
could also be loaded in the main matlab file and it would 
ultimately give us the SNR that allowed us to estimate 
which prediction order or bits per sample are appropriate to 
re-generate the signal as close to the original signal as 
possible. 

SNR was calculated by the given formula : 



Reconstructed signal after Dequantization 



J I I I L 


4. Experimental Results 

In the experimental section I basically wanted to choose a 
suitable prediction order and the number of bits/sample to 
make a good compromise between the compression and the 
sound quality. We did that by comparing the SNR values 
of obtained from the following equation. 

Hy[n] 2 

SNR = 10log 10 \ (6) 

Ziyln] - r[n\) 2 


Which is just a modified version of the following equation: 

SNR = -10 log 10 ^^ (7) 


I started the procedure by keeping K-l and calculated the 
values of both SNR_r and SNR_e for 1< N < 5. 


Table 3: SNR Comparison for k=l and 1<N<5 


K 

N 

SNR r 

1 

1 

-22.1855 

1 

2 

-16.1649 

1 

3 

-10.2901 

1 

4 

-6.1562 

1 

5 

-4.6608 

1 

6 

-4.3322 


I observed from the table that as I increased the prediction 
order N, the SNR_r increased, which is desirable as the 
system was able to produce a sound quality much similar to 
the original signal. The system was thought to have a better 
Compression i.e. a better prediction then. As I saw that the 
system was getting a better sound quality as well as a better 
compression for both the increasing K and N I constructed 
the following table to get the optimal SNR_r and SNR_e 
values at the range of 1 < K < 11 and 1< N < 10. 

2 

1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 

0 1 2 3 4 5 6 

xIO 4 

Figure 7: Transmitted Signal for K = 4 N = 4 


transmitted signal 



Table 4: SNR Comparison for 1<K<1 1 and 1<N<1 1 
(laughter) 


K 

N 

SNR r 

1 

1 

-22.1855 

2 

2 

-11.6936 

3 

3 

-6.6188 

4 

4 

-2.9897 

5 

5 

-2.2874 

6 

6 

-1.9291 

7 

7 

-1.7244 

8 

8 

-1.6341 

9 

9 

-1.6012 

10 

10 

-1.4655 

11 

11 

-1.4442 
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transmitted signal 


5 

4.5 



0 1 2 3 4 5 6 

xIO 4 

Figure 9: Transmitted Signal for K = 5 N = 5 



-i 1 1 1 u 1 

0 1 2 3 4 5 6 

xIO 4 

Figure 10: Original Vs. Decoded signal for K = 5 N = 5 


transmitted signal 




0 1 2 3 4 5 6 


xIO 4 

Figure 12: Original Vs. Decoded signal for K = 6 N = 6 


transmitted signal 
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Figure 14: Original Vs. Decoded signal for K = 7 N = 7 


Table 5: SNR Comparison for 1<K<10 and 1<N<10 (Handel) 


K 

N 

SNR r 

1 

1 

-24.4608 

2 

2 

-14.7124 

3 

3 

-8.6330 

4 

4 

-4.6088 

5 

5 

-4.3354 

6 

6 

-3.4003 

7 

7 

-2.9174 

8 

8 

-2.3727 

9 

9 

-2.1369 

10 

10 

-2.0144 


Table 6: SNR Comparison for 1<K<10 and 1<N<10 (Train) 


K 

N 

SNR r 

1 

1 

-19.1465 

2 

2 

-8.2138 

3 

3 

-3.2407 

4 

4 

-1.3343 

5 

5 

-1.3125 

6 

6 

-1.4291 

7 

7 

-1.5104 

8 

8 

-1.0929 

9 

9 

-1.1158 

10 

10 

-0.7450 



Quantization Level (1<K<10) and Prediction Order (1<N<10) 

Figure 15: Signal to Noise ratios for both r and e against ( 1 <k< 1 1) and 
(1<N<1 1) for laughter signal 



Quantization Level (1 <K<1 0) and Prediction Order (1<N<10) 

Figure 16: Signal to Noise ratios for both r and e against (l<k<10) and 
(1<N<10) for handel signal 



Quantization Level (1<K<10) and Prediction Order (1<N<10) 

Figure 17: Signal to Noise ratios for both r and e against (l<k<10) and 
(1<N<10) for train signal 

From the graph above it can be said that SNR r is maximum 
for both K = 10 and as well as N = 10. Whereas the 
prediction error is lowest for lowest for K = 10 and N = 10. 
So it can be said that when the signal is quantized using 
increasing number of bits per sample as well as a higher 
prediction order it can have a proper audible reconstruction 
of the signal as well as a better compression due to 
decreasing e[n]. 



6 



Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 16, Issue 2, ICGST, Delaware, USA, December 2016 


It can also be seen that both the SNRs approximately 
intersect each other in a point when both (K, N) = (4, 4). So 
the optimum point is seen somewhere around (K, N) = (4, 
4) for both the laughter and the handel signals. But for train 
signal the optimum point will be at (K, N) = (5, 5). 


Table 7: Optimum SNR for trade-off between K(2-l 1) and 
N(l-6) for TRAIN Signal 


N 

K 

0 

1 

2 

3 

4 

5 

6 

2 

-19.59 

-13.71 

-8.21 

-3.82 

-1.84 

-1.61 

-1.61 

3 

-18.56 

-12.72 

-7.35 

-3.24 

-1.63 

-1.44 

-1.41 

4 

-17.40 

-11.60 

-6.41 

-2.66 

-1.33 

-1.32 

-1.32 

5 

-17.41 

-11.62 

-6.42 

-2.67 

-2.67 

-1.31 

-1.32 

6 

-16.10 

-10.61 

-5.42 

-2.10 

-1.20 

-1.44 

-1.43 

7 

-13.42 

-7.96 

-3.64 

-1.24 

-0.70 

-1.38 

-1.49 

8 

-9.51 

-4.75 

-1.76 

-0.51 

-0.25 

-1.02 

-1.01 

9 

-8.16 

-3.78 

-1.30 

-0.37 

-0.12 

-0.81 

-1.02 

10 

-5.80 

-2.31 

-0.70 

-0.19 

-0.04 

-0.75 

-0.68 

11 

-11.56 

-6.37 

-2.64 

-0.83 

-0.22 

-0.65 

-0.78 


Table 8: Optimum SNR for trade-off between K(2-l 1) and 
N(7-10) for Train Signal 


N 

K 

7 

8 

9 

10 

2 

-1.62 

-1.63 

-1.63 

-1.63 

3 

-1.42 

-1.43 

-1.43 

-1.42 

4 

-1.33 

-1.33 

-1.33 

-1.33 

5 

-1.32 

-1.32 

-1.33 

-1.33 

6 

-1.43 

-1.43 

-1.42 

-1.42 

7 

-1.51 

-1.59 

-1.53 

-1.54 

8 

-1.05 

-1.09 

-1.09 

-1.09 

9 

-1.03 

-1.10 

-1.11 

-1.12 

10 

-0.73 

-0.75 

-0.75 

-0.74 

11 

-0.79 

-0.72 

-0.76 

-0.77 


So it can be seen from the table above that the optimum at 
Prediction Order, N = 4 and Uniform quantization of Bits 
per sample, K = 10 we got the highest SNR of -0.04 for train 
signal. Similarly the optimum SNR for trade-off between K 
and N for laughter as well as handel signals could be found. 
As these optimal values of K and N are determined for the 
three signals train, laughter and handel. Similarly this 
custom built GUI and MATLAB functions could be applied 
to any functions with after loading them as raw data inside 
MATLAB. So as discussed earlier in the introduction 
section the applications mentioned in the papers [2, 3, 4, 5, 
6] could use the simulation process of this paper and 
increase the efficiency of their applications by using the 
optimal values of the bits per sample of the uniform 
quantization process in their applications. 


5. Conclusion 

This system successfully recreated both the encoder and 
the decoder of open-loop DPCM by applying some custom 
built MATLAB functions. It was basically used to find the 
trade-offs between the bits per sample of the uniform 
quantization process and the prediction order by comparing 
their SNRs. The results that were obtained could be used in 
variety of applications such as compression, de-noising, 
video coding, prediction, interpolation, and 
communications. The proposed technique of modeling and 
simulation of this paper is restricted to the application to 
DPCM simulation only. For an improvement to this paper, 
the proposed modeling and simulation technique could be 
applied to wide variety of digital filters like Kalman filter, 
wiener filter, LMS, RLS, AR, ARX, ARMA, ARMAX, HR, 
FIR, FFT, IFT, DCT, DST etc. 
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Abstract 

This paper develops a multilevel structure based on the 
modeling and control of a Photovoltaic (PV) application 
connected to a three-phase grid. The DC-DC boost 
converter boosts the DC voltage of the PV array to the 
Maximum Power Point (MPP) using a Maximum Power 
Point Tracking (MPPT) based on the incremental 
conductance method. For synchronization, a three-phase 
inverter is necessary to connect the PV system to the grid. 
The control system modeling of three phase inverter was 
carried out under Matlab / Simulink environment using 
the proportional integral controllers, pulse width 
modulation (PWM) three-level pulse generator technique 
and park transformation. The Phase Locked Loop (PLL) 
is used to synchronize the inverter current angle with the 
angle of the grid voltage to reach a power factor as close 
to one as possible. As a final point, different types of 
fault condition are created, simulated and analyzed on the 
grid side. 

Keywords: photovoltaic (PV), Maximum power point 
tracking (MPPT), pulse width modulation (PWM), 
voltage source converter (VSC), fault analysis 

Nomenclature 

I D The diode current (A) 

V d The diode voltage (V ) 

I sat The diode saturation current ( A) 

y ^ The thermal voltage 

T The cell temperature ( K ) 

K The Boltzmann constant, equal to 

1.3806 c -23 (•/•£'’) 

q The electron charge, equal to 


1.6022e~ 19 (C) 

Q d The diode quality factor 

j? The shunt resistance 

1X sh 

j? The series resistance 

At The number of series-connected cells 

1 V cell j , 

per module 

yy The number of series-connected 

module per string 

yy The number of parallel strings 

pp 

&base The base v °l ta g e ( ^rms line-to-line). 

p The inductive three-phase short circuit 

power in (VA) 

/ The frequency in (//Z) 

X The reactance in (Q) 

(X/Fty The ratio 

1. Introduction 

Today, the energy demand in the world is in full 
evolution. However, the stock of oil in our planet will 
soon be exhausted. The global warming of the planet 
today becomes more serious due to the greenhouse effect. 
Some emissions of these gases are due to human activity. 
To reduce these emissions and ensure energy security, 
renewable energy is enhanced in the production of 
electrical energy while preserving our fossil deposits (oil 
and gas) as long as possible with a lesser pollution and 
controlled. 

In recent years, there has been a growing interest in PV 
solar energy as it is a potential source for electricity 
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generation with minimal environmental impact [1]. This 
energy is abundant, distributed through the earth, 
pollution free and clean and noise-free source of energy 

[2] -[3], 

Direct exploitation of solar energy falls within three 
distinct technologies: one produces calories, it is solar 
thermal, the second is thermoelectric generators, and the 
third produces electricity, it is photovoltaic (PV) solar 
energy, which is the subject of this paper. The growing 
interest in photovoltaic systems requires the growth of 
the activities of research and development in various 
aspects such as Maximum Power Point Tracking (MPPT) 
[4], PV arrays, anti-Islanding protection, stability and 
reliability, the electronic interface of power quality and 
power. The MPPT is used to control variations on the 
curve current voltage of the cells. Thus the MPPT system 
must track the point of optimal functioning of the PV 
arrays. The purpose of this technique is to ensure that the 
system operates very close to the optimal operating point. 
In regards to power electronic converter interface PV 
array to the grid, a DC-DC boost converter and a three- 
phase three-level Voltage Source Converter (VSC) are a 
topology widely used to this day. We used 
Matlab/Simulink environment to model the control 
system of three phase inverter using the proportional 
integral controllers, pulse width modulation (PWM) 
three-level pulse generator technique and park 
transformation. To get the power factor as close as 
possible to one, we use the Phase Locked Loop to 
synchronize the inverter current angle with the angle of 
the grid voltage. 

Finally, fault analyses is carried by creating a line-to- 
ground (LG), line-to-line (LL), double line-to-ground 
(LLG) fault in the network and observe the ac voltage, ac 
current and ac power waveforms at the grid side with and 
without different types of fault condition. The remainder 
of the paper is organized as follows: Section (2) focuses 
on description of the grid connected PV system. Section 

(3) emphasizes on PV array model. Section (4) is 
dedicated to the control system of PV solar systems. 
Section (5) presents the grid model. Section (6) provides 
a calculation of asymmetrical short-circuit current: 
Method of symmetrical components. Our results are 
analyzed and discussed in section (7). The main 
conclusions are drawn in section (8). 

2. System description 

The principal schematic diagram of the grid connected 
PV system is presented in figure l.The PV module is 
usually composed of a number of solar cells with 
identical characteristics. Similarly, a number of PV 
modules are assembled and interconnected to build a PV 
array. The PV array generate 72.6 KW under standard 

test conditions ^25° C, 1000W / m 2 ^ . The voltage 

level of the PV array 243V is boosted to 400 V via a 
DC-DC boost converter using MPPT based on the 
incremental conductance method. This dc power is 
converted into ac power using a three level bridge 
inverter. This inverter has three legs with two switches in 
each leg. The switching is performed by carrier-based or 


space-vector based pulse width modulation. The 
converter produces an output waveform very close to a 
sinusoidal wave with extremely low harmonic distortion 
levels. The output Zy filter is connected to remove high 

switching frequency components from output current of 
inverter [5]. A delta-star step up transformer configuration 
is generally used in grid connected system because the 
third harmonic will get circulated in delta and does not 
enter in the grid. To produce the pulses for the three-level 
bridge inverter, a dq control uses a reference frame 
transformation abc to dq frame which transforms the 
grid current and grid voltage PLL into dq frame. 

3. PV array model 

The basic equation that describes the I-V characteristics 
of the PV model is given by the following equation [2]- 
[ 6 ]: 

I D =I sal [cx P (V d /V T )-l] (l) 



where: 

p, 

q 

4. Control system 

The system has multiple control blocks that work 
together to ensure maximum power extraction from the 
PV array and then converting it to AC power to be 
injected into the grid. In this section, we discuss the 
details of each control block in the system. 

4.1. Control of DC-DC boost converter 

The DC-DC boost converter boosts the DC voltage of the 
PV array to the MPP using an MPPT based on the 
incremental conductance method. This method computes 
the maximum power point by comparison of the 

incremental conductance (A/ / AV) to the array 
conductance (//V) . When these two are the 

same ( I IV = A7 / A V) , the output voltage is the MPP 

voltage. The controller maintains this voltage until the 
irradiation changes and the process is repeated. The 
incremental conductance is a technique that takes 
advantage of the fact that the slope of the power- voltage 
curve is zero at the maximum power point. The slope of 
the power- voltage curve is positive at the left of the MPP 
and negative at the right of the MPP [7] as given by: 
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Figure 1. Schematic diagram of the grid connected PV system 


dP / dV = 0, at MPP 
<dPldV> 0, left of MPP 
dP / dV < 0, right of MPP 

since 

dP d(lV ) dl A I 

dV dV dV AV 

Therefore the equation (4) can be rewritten as: 

A//AV = -I/V, at MPP 
<AI/AV>-I/V, left of MPP (6) 

AI/AV<-I/V, right of MPP 

We can find the MPP by comparing the instantaneous 
conductance (//V) to the incremental conductance 

(A 77 AV ) as shown in figure 2. This idea is presented 
in [8]. 

This technique has an advantage over the perturbation 
and observation (P&O) method because it can 

determine when you reach the MPP without having to 
oscillate around this value. It can also perform MPPT 
under rapidly increasing and decreasing irradiance 
conditions with higher accuracy than the perturbation and 

observation (P&O) method. 



Figure 2. Flowchart for an Incremental Conductance Tracking System 


4.2. Control of voltage source converter (VSC) 

The three level VSC regulates DC bus voltage and keeps 
unity power factor. The control system uses two control 
loops: an external control loop which regulates DC link 
voltage and an internal control loop which regulates I d 

and / grid currents (active and reactive current 
components). I d current reference is the output of the 
DC voltage external controller. / current reference is 
set to zero in order to maintain unity power factor. V d 
and V voltage outputs of the current controller are 
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converted to three modulating signals U ref abc used by 

the PWM three level pulse generator [9]. The 
determination of the grid voltage vector angle is done by 
a phase-locked loop (PLL) algorithm. 

The grid voltage vector is used in a feed forward 
loop to compensate for the grid harmonics [10]. The 
voltage and power equations by the inverter in 
d — q synchronous reference frame are represented 
below: 

v: = RI. +V, - coLI + L— 1, 

d d d “ dt d 

(7) 

V; = RI q+ V q -vLI d+ L-I q 

P = v d l d +V q l q 

( 8 ) 

Q = ~vj d +v q i d 


been developed and presented in many papers up to now 
[15]- [16]. The most common method for three-phase 
system is that based on the representation in the 
referential Park of the measured voltages on the network. 
The principle is to consider the equations of the three- 
phase network in a reference Park whose orientation is 
such that V q =0 [17]. 

The phase locked loop achieves the slaving an arbitrary 
phase angle to the reference angle. The simulink model 
for a three-phase PLL is shown in figure 4. 

As it can be observed, a PID controller is used in 
order to reduce the error between the reference and 
measured values of V . Deduce dq components from 

abc signals by performing an abc to a (3 Clarcke 
transformation in a fixed reference frame, then perform 
an a/3 to dq transformation in a rotating reference 

frame to obtain the measured value V . 


The voltages presented by equation (7) in the 
d — q reference are coupled. Figure 3 shows how 
compensation is made for the coupling of the axes d and 
q. From the voltages equation by the inverter in d-q 
synchronous reference frame, we obtain the following 
decoupled voltages equation [11]: 

W; = RI d + V d -coLI q + K p (/* -I d ) + Kl \(f d ~h)dt {9) 

\v;=RI q+ V q+ wLI d+ K p {ll-I q ) + K,\{ll-I q )dt 
The bloc diagram is presented in figure3. 



Figure 3. Compensation of the cross coupling between axes 


4.3. Filter inductance L j 


To reduce the switching losses in the power components 
of the inverter, and aimed at improving the load current 
(injected current), was added an inductance [12]. The 
value of L f is design based on current ripple. Typically 


the ripple current can be chosen as 15% - 25% of rated 
current. The maximum current ripple [13] -[14] can be 
resulting as following: 


Air 


1 V, 


dc 


8 L f f s 


= 20 % • 


( 10 ) 


'fJsw 


5. Grid model 

The grid was modeled as a Y connected a balanced three- 
phase voltage source with an internal R g - L g impedance 
[18]. You can specify the source internal resistance and 
inductance either directly by entering R g and L g values 

or indirectly by specifying the source inductive short- 
circuit level and X / R ratio. The internal inductance 
and resistance are defined as: 


P sc 2 ^f 

(X/R) (X/R) 


( 11 ) 

( 12 ) 


6. Calculation of asymmetrical short-circuit 
currents: Method of symmetrical 
components 

Any type of sinusoidal waveform (voltage, current, 
flux ...) can be written in terms of its symmetrical 
components as follows [19]: 


X' 


i 

>< 
i 


* 

ii 

x 2 

X. 


i 

X 

0 

1 


with 

f i 

1 

ij 

.In 

[ p ] = 

a 2 

a 

i 

and a- e 3 


v a 

a 2 

b 



(13) 


4.4. Phase locked loop (PLL) 

The PLL command is necessary in the control of the 
converters. Different methods to extract phase angle have 
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Figure 4. The Simulink model for a three-phase PLL 


The sequences positive, negative and zero are determined 
by: 
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7. Simulation results 

In this section we present the simulation results of the 
proposed grid connected PV system presented in figure 1 . 

7.1. Simulation ofPV system 

The PV array consists of 66 strings of 5 series-connected 
modules connected in parallel. The total power 

is (66 * 5 * 220W = 12.6KW ) . Manufacturer 

specifications for one module are shown in Appendix 
(Table 1). 

Figures 5, 6, 7 and 8 shows the output 
characteristics of PV arrays. These curves are nonlinear 
and are crucially influenced by solar radiation and 
temperature. 



Figure 5. V-I characteristics of array under constant temperature and 
different irradiance 



Figure 6. V-P characteristics of array under constant temperature and 
different irradiance 




— T=25°C 
— T=35 8 C 
T=45°C 
T=55°C 
— T=65 # C 


150 200 250 300 350 

Voltage (V) 


At constant irradiance G=100t) \\7nr 


150 200 

Voltage (>) 

Figure 7. V-I characteristics of array under constant irradiance and 
different temperature 

n x10l 


T=25 C 


At constant irradiance G=1000 \\ /nr 


T=35°C 


T=45°C 


T=55 : C 


T=65 C 


Figure 8. V-P characteristics of array under constant irradiance and 
different temperature 
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7.2. Simulation of grid connected PV system without fault 

The grid connected PV system without fault in figure 1 is 
simulated by using MATLAB/Simulink software 
environment in order to validate the control methodology 
discussed earlier. Parameters of the test system are 
summarized in Appendix (Table 2). 

For all simulations, the temperature is assumed constant 

at 25 °C . Under steady-state condition, (i.e. constant 
solar irradiance G = 100CHV/m 2 ), some simulation 
results have been chosen to illustrate the system response 
at standard test condition (STC). 


15 =i 

-la 

—15 



3 * W* 

Vi V6 Vl 



3 o“ 0.02 0.04 0.M COS 0.1 

Time (s] 

Figure 9 (b). The three phase grid voltage waveform 


a 

i 3 





Figure 9 (e). Inverter output voltage 

700 




Figure 9 (g). Dynamic response of synchronization of the injected 
current 
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Figure 9 (h). The total harmonic current distortion in frequency domain 


The three phase grid current and grid voltage waveforms 
are exposed in Figures 9(a) & 9(b) respectively. As can 
be seen, the system response under STC is very good 
where the steady-state conditions are reached in few AC 
periods. The dq components of the injected current are 
shown in Figure 9(c) which follow closely the references 
l d and/^ . The real power and reactive power injected 

from the PV system are shown in Figure 9(d). The 
injected real power was around 1\KW while reactive 
power was set to zero. To eliminate the high frequency 
harmonic in Figure 9(e), low-pass filter is connected. The 
DC link voltage was controlled at a steady state value of 
about 400V as shown in Figure 9(f). The injected 
currents must be synchronized with the grid voltages. So, 
the PLL has been used to satisfy this condition. The 
dynamic response of synchronization of the injected 
current is shown in Figure 9(g). The total harmonic 
distortion (THD) of the current in frequency domain is 
shown in Figure 9(h). The THD of current injected into 
the grid is 2.47% which satisfies IEEE 519 standard 
limits. 

7.3. Effect of irradiation 

In order to investigate the effectiveness of the control 
algorithms, a dynamic simulation is carried out with solar 
irradiance suddenly dropped from lOOOFF / m 2 to 
2501V / m 2 at 0.5 sec as shown in Figure 10(a). 

The effect of irradiation on PV power generation is 
demonstrated in figure 10. The PV array output current, 
voltage and power during simulation run are shown in 
Figures 10(b), 10(c) and 10(d) respectively. It may be 
noted that the effect of solar irradiance on the array 
terminal voltage is small compared to its effect on output 
current. 

The simulation results of irradiation effect on PV 
connected to grid are illustrated in figure 1 1 . The active 
and reactive powers are shown in Figure 11(a). The 
reactive power drops from 71 KW to 14 KW after 
0.5 sec when the reactive power is maintained equal to 
zero. From the simulation curve in figure 11(b), can be 
distinguished that the grid current amplitude decreases in 
response due to the sudden change of the PV generated 
power. 
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Figure 10 (a). The solar irradiance 
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Figure 10 (d). The PV array power 



15 



Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 16, Issue 2, ICGST, Delaware, USA, December 2016 



Figure 1 1 (a). Active and reactive power injected to the grid 



0.1 0.2 0.3 0.4 0.5 0.6 07 0.8 

Time- (i| 

Figure 1 1 (b). The grid current 


7.4. Simulation of grid-connected PV system with fault 

The system is simulated at temperature T = 25 °C and 
irradiance G = 1000VL / m 2 . We applying a fault on the 
grid side at £ = 0.55sec and the fault duration 
is 0.05 sec . 

• Line-to- ground (LG) fault 


„xia* 


V4 Vb ve 
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Time 


Figure 12. The grid voltage with LG fault 
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Figure 13. The grid current with LG fault 


e 100 



Due to LG fault, the voltage V a drops to zero as shown 

in figure 12. Application of fault on the grid side has 
resulted in oscillations in current. However, oscillations 
are damped as soon as the fault is cleared, and the system 
tracks the reference current in less than 0.2 sec. The 
injected current into the grid is no longer sinusoidal. This 
disconnection creates an increase over rated currents of 
the three phases currents simultaneously as shown in 
figure 13, this case study clearly shows an increase equal 
to 16.16 times of the rated current for phase A, whereas it 
is equal to 17.44 times for phase B and 14.24 times for 
phase C (see Table 3 in Appendix). In addition, this 
voltage failure causes a slight decrease in the level of 
active and reactive powers. In figure 14, the decrease of 
the active power is about 0.93 times of the rated power 
(see Table 3 in Appendix). 


Figure 14. The active and reactive power with LG fault 

• Line-to- line (LL) fault 

When applying the LL fault, it can be observed from 
figure 15 that two phases A and B are in short-circuit 
during a period of fault. As a result, application of fault 
on the grid side has resulted in oscillations in current so 
that an enhances observed as shown in figure 16; this 
injected current into the grid is no longer sinusoidal and 
this boost is less important than the LG fault. For phase A, 
an increase of 2.17 times is reached, while for phase B an 
increase of 1.81 times and for phase C an increase of 2.5 
times are obtained (see Table 3 in Appendix). This case 
study clearly shows an decrease in active power is more 
important for this type of fault than for the LG fault as 
shown in figure 17, and this decrease is 0.86 times (see 
Table 3 in Appendix). 
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Figure 18. The grid voltage with LLG fault 


Figure 19. The grid current with LLG fault 
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Figure 15. The grid voltage with LL fault 
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Figure 16. The grid current with LL fault 
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Figure 17. The active and reactive power with LL fault 


• Line-to- line to ground (LLG) fault 

The different simulations of LLG fault are shown in 
the following figures. As shown in figure 18, phases 
A and B are grounded. This fault applied in the grid 
side causes an increase of the phase currents. 
However, this increase is damped as soon as the fault 
is cleared and the system tracks the reference current. 
This increase is allocated as follows: phase A- 10.72 
times, phase B- 14.84 times and Phase C- 12.62 
times (see Table 3 in Appendix). This increase is 
smaller than the LG fault but at the same time larger 
than the LL fault as shown in figure 19. During fault 
period, the active and reactive power decrease in 
figure 20 is less than LG and LL faults. The decrease 
in active power is 0.71 times (see Table 3 in 
Appendix). 


£150- 



Time (a) 


Figure 20. The active and reactive power with LLG fault 

8. Conclusion 

In this paper, a grid-connected PV system based on 
MATLAB has been proposed. The main objectives were 
to achieve maximum power output from the PV array and 
to inject a high quality AC current into the grid. 

Based on the modeling system, MPPT, steady state and 
dynamic performances, and fault analysis are studied. 
The PLL has been designed for grid synchronization and 
it effectively synchronizes the inverter voltage and 
frequency with the grid voltage and frequency. The 
steady state and dynamic performances of the proposed 
system have been tested under constant solar irradiation 
and fast changing solar irradiation. The presented 
simulation results show a good system behavior under all 
these case studies. 

Fault analysis on grid side have been performed for 
various fault conditions like LG, LL, LLG faults. In case 
of fault, to become stable at nominal frequency, it is 
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observed that the system tracks the reference current in 
less than 0.2 sec. Application of fault on the grid side has 
resulted in current oscillations. The injected current into 
the grid is no longer sinusoidal. Disconnections create an 
increase over rated currents of the three phase currents 
simultaneously. In addition, this voltage failure causes a 
slight decrease in the level of active and reactive powers. 
The results of simulations show that the interruption due 
to phase to ground faults are noticeably more dangerous 
than the line to line faults. The decreasing of active 
power depends on the number of phases in fault. Nearly 
all PV systems will soon be required to have a fault 
protective device that will minimize the damages in PV 
arrays. 

9. Appendix 

• Manufacturer specifications for one module are 
shown in table 1 


TABLE 1 

ELECTRICAL CHARACTERISTIC FOR ONE 


MODULE 

Symbol 

Manufacturer 

Rating 

values 

P 

max 

Maximum power 

220W 

N 

1 V cells 

Number of series-connected 
cells 

96 

Voc 

Open-circuit voltage 

59.2618V 

I sc 

Short-circuit current 

5.09261A 

V 

mp 

Voltage at maximum power 

48.3159V 

1 mp 

Current at maximum power 

4.54758A 


• System parameters used in simulation are shown in 
table 2 


TABLE 2 

SYSTEM PARAMETERS USED IN SIMULATION 


Components 

Rating value 

DC link: 

Capacitance 

\2000 pF 

Reference V dc 

400V 

Passive L Filter: 

250//// 

Inductance: L 

Resonance frequency 

50HZ 

R (in series with L) 

0.0020 


Step-up 

Isolation Transformer: 

100/TVA 

Apparent power: S 

Primary side 

25 KV ( L-Lrms ) 

Secondary side 

230V ( L-Lrms ) 

3 phase AC distribution grid: 

25 KV 

Phase-to-phase rms voltage 

Frequency 

50 HZ 

3 phase short circuit level 

100MVA 


• Phase currents and active power variations during 
different fault conditions are shown in table 3 

TABLE 3 

PHASE CURRENTS AND ACTIVE POWER 
VARIATIONS DURING DIFFERENT FAULT 
CONDITIONS 


FAULT 

PHASE 

PHASE 

PHASE 

ACTIVE 

A 

B 

C 

POWER 

LG 

1616% 

1744 % 

1424 % 

93% 

LL 

217 % 

181.4% 

250 % 

86% 

LLG 

1072 % 

1484 % 

1262 % 

71 
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Abstract 

Sensor less Control of PMSM with FOC using MRAS for 
speed estimation is implemented through Popov’s Hyper 
stability criterion. The output of the adaptation mechanism 
MRAS is the estimated speed quantity using PI controller. 
The aim of the proposed sensorless control is to improve 
performance and robustness of PMSM drive under load 
and speed variations. The PMSM drive is tested for three 
different cases; with balance three phase supply, FOC and 
with MRAS using PI controller. The effectiveness and 
validity of the proposed control approach is verified by 
simulation results through MATLAB/Simulink/Sim 
Power System environment. It is found that the 
performance of the PMSM drive is better and also the 
torque, flux ripples are quite less with MRAS. 

Keywords Permanent magnet synchronous motor 
(PMSM), Field oriented control (FOC), Model reference 
adaptive system (MRAS). 


Nomenclature 

V ds direct axis stator voltage, V 

V qs quadrature axis stator voltage, V 

ids direct axis stator current, A 

i qs quadrature axis stator current, A 

R s Stator armature resistance, Q 

p no. of poles 

Lds direct axis inductances, H 

L qs quadrature axis inductances, H 
Ads direct axis flux linkage, wb 
A qs quadrature axis flux linkage, wb 

Af magnetic flux linkage, wb 

C0e rotor speed in electrical, rpm 

C0 m mechanical speed, rad/s 

J moment of inertia, kg.m 2 

B Viscous Friction Co-efficient, Nm/rad/s 

P Differential operator 

T e electromagnetic torque, Nm 

T l load torque, Nm 


1. Introduction 

As PMSM Motor is made up of rare earth and neodymium 
boron magnets, it has been widely used in high 
performance variable speed industrial applications. In 
this motor, Permanent Magnets are placed in the rotor; 
because of absence of windings in the rotor; rotor copper 
losses are zero. Due to these advantages, this motor offers 
high effiency, high torque to inertia, high Power density 
[1]. PMSM are preferred in applications where it requires 
fast torque response and high performance operations [2] 
such as robotics, electrical vehicles, and servo applications 
[3, 4]. Both Permanent Magnet Synchronous Motor 
(PMSM) and Brushless DC motor (BLDCM) have 
permanent magnets on the rotor and require alternating 
stator currents to produce constant torque. But the torque 
ripples of the BLDCM are higher than that of the PMSM 
[ 1 ]. 

The scalar control is an open loop speed control method, 
doesn’t provide a possibility to control the currents during 
various operation cycles like the start-up and loading 
conditions. The second method, vector control ( FOC & 
DTC) is the most common method of speed control 
method for AC drives due to their dynamic response [5]. 
Comparison between FOC and DTC are based on various 
criteria including basic control characteristics, dynamic 
performance, and implementation complexity. The 
simulation and evaluation of both control strategies are 
performed using actual parameters of Permanent Magnet 
Synchronous Motor fed by an IGBT PWM inverter [6]. It 
presents a method of estimating simultaneously the motor 
speed and the rotor resistance of an induction motor by 
superimposing ac components on the field current 
command. The validity of the proposed method is verified 
by simulation and experimentation [7] . 

The structure and the control methods of PMSM are 
analyzed and simulation is realized using conventional 
direct torque control (DTC) method. It can be applied for 
PMSM and is reliable in a wide speed range applications 
[8]. Dynamic permanent magnet Flux Estimation of 
Permanent Magnet Synchronous Machines explains 
clearly about the rotor flux estimation in permanent 
magnet synchronous motor [9]. A Direct Torque 
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Controlled Interior Permanent Magnet Synchronous 
Motor Drive without a Speed Sensor has been presented 
with a speed control system for an IPM motor with an 
inner DTC which does not require a mechanical position 
sensor [10]. 

The problem of efficiency optimization in vector- 
controlled interior permanent magnet synchronous motor 
is investigated. It also describes a method for minimizing 
the losses in vector controlled interior permanent magnet 
synchronous motor drives [1 1]. A new control scheme for 
wide speed range operation of interior permanent magnet 
synchronous motor drives, where both torque and stator 
flux linkage are directly controlled is proposed in [12] .It is 
proved that the dynamic performance of permanent 
magnet synchronous motor with vector control is better 
[13]. But the main advantage of DTC is that there is no 
requirement for transformation, and current regulator, but 
it gives high torque ripples [14]. The basic principle of 
DTC is to directly select stator voltage vectors according 
to difference between the actual & reference torque and 
stator flux linkages. In this paper along with stator 
resistance estimator is used [15]. In the DTC torque and 
flux ripples are high; to overcome this problem fuzzy 
controller is used. Using hysteresis band the torque ripples 
are reduced and MRAS technique is used for estimate the 
motor speed [16]. To trace the error between reference and 
actual values, the conventional controller PI, PID is 
implemented. To have better performance and accuracy in 
outputs, Model Reference Adaptive Control system 
MRAS is adapted [17]. The performance of the PMSM 
can further improved by implementing SVM control 
strategy along with DTC and Multi Level Inverters using 
Fuzzy Logic Controllers [18] [19] [20] [21] [22]. The 
hardware implementation of PMSM with PI Controller 
using FPGA is discussed in [23] using two level inverter 
with SVM control. 

In this paper, Sensor less Control of PMSM with FOC 
using MRAS is implemented. The rotor speed is estimated 
using MRAS with PI controller. The paper is organized as 
follows: Section II revolves around the Field Oriented 
Control of Permanent Magnet Synchronous Motor. 
Section III explains the implementation of FOC for PMSM 
drive with the proposed MRAS. Section IV focuses on the 
simulation results for the dynamic performance of PMSM. 
Section V summarises the conclusions drawn from the 
work with future recommendations. 


2. Field Oriented Control of Permanent 
Magnet Synchronous Motor 

The principle idea of FOC is to introduce decoupling 
between field and torque producing component. This 
makes PMSM to behave like a separately excited DC 
motor. If the magnets are placed inside the rotor then L q > 
L d , otherwise if the magnets are placed on the surface of 
a rotor then L q = L d . 

The modeling of PMSM motor in rotor reference frame is 
given from Eqns (1) to (8). The proposed block diagram 
of the work is given in Figure 1 . 


^qs — Rs^qs T P^qs T ®r^ds 

(1) 

Vds — R s ids T P^-ds — ®r^-qs 

(2) 

^ds “ LdsUs + 

(3) 

^-qs — Lqs^qs 

(4) 

v ds — P^sUs T Pds^ids — Lqs^riqs 

(5) 

v qs — Ps^qs T Lqs ~[[iqs Pqs^rUs 

(6) 

Te — ^K^ds — L , q S )i d siqs ^-fiqs] 

(7) 

J^f+Bco m + T L = T e 

(8) 


FOC 



3. Implementation of FOC for PMSM drive 
with the proposed MRAS using PI 
Controller 

The Field Oriented Controlled electrical drives have rapid 
expansion in recent years, due to the achievements 
obtained in semiconductors for both power and signal 
electronics. This paper proposes to describe Field oriented 
control (FOC) of Permanent Magnet Synchronous motor 
(PMSM) using MRAS speed observer for sensor less 
control of drive. 

This paper proposes a sensor less speed control based on 
MRAS, which is based on the comparison between outputs 
of two estimators. This estimated error is used to drive 
PMSM with a suitable adaptive mechanism to estimate the 
error speed. In this work MRAS is based on i d i q of a 
PMSM, since FOC is used as a speed controller. 

The general block diagram of MRAS is given in Figure 2. 
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i — 1 



Figure 2.Block Diagram of MRAS 

The main objective of MRAS is to estimate the rotor speed. 
Hence first i d i q are represented in a state variable form is 
given by Eqn (9) 



According to Eqn (9), the following Eqn (10) is derived 


id+ t 


. iq . 



R s Lq' 

C0 P — 

L d L d 

Ld R s 

-co P — 

Lq Lq _ 


id+ ^ 


+ 


V d I R S y 


( 10 ) 

From Eqn (10) the reference currents and voltages are 
obtained as given in Eqn (11) 

id = id ,iq = i q ,v5 = Vd+^Xf.V* =v q (11) 

Using Eqn (11), Eqn (10) can be converted to Eqn (12) 


4 


-5 * co p Ll' 


L d 
L d 
u e i 


' L d 


- r 

n e. n 


L d U 

■S’!. 


( 12 ) 


According to Eqn (12), the state equation of adjustable 
model of PMSM with speed angle as the adjustable 
parameter is represented in Eqn(13): 



(13) 


For SPM, Ld=Lq=Ls, so the adaptive speed mechanism can 
be simplified as follow: 

®e = ( Kp + j) [idlq - iqld “ 0q “ Q] + ^(0) 

(14) 

Block Diagram of FOC of PMSM with MRAS is given in 
Figure 3. 
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Figure 3. Block Diagram of FOC of PMSM with MRAS 

4. Design of PI Controller - Conventional 
Controller 

For process control in industrial automation the 
conventional PI controllers are extensively used because 
of its simple process, design, maintenance, economical, 
effectiveness, less number of parameter tuning and good 
performance for different operating conditions. Also it 
does not require the complete knowledge of the system. 
Hence these controllers are mostly preferred in industrial 
applications. [24] [25] [26]. The integral term in a PI 
controller causes the steady-state error to reduce to zero. 

Speed controller plays an significant role in achieving 
aspiring dynamic characteristics of AC drives (PMSM). In 
this work, PI controller is considered in speed control loop 
[27] to investigate the dynamic response of a drive for 
different set speeds and load conditions. This is achievable 
by tuning the parameters of a PI controller, K p and Ki [25]. 

The stability of a IPM and SPM drive system are 
analyzed as a linearized error model using PI controller. 
The general block diagram of the PI speed controller is 
shown in Figure 4.. 



Figure 4 Block diagram of PI speed controller 
The output of a speed controller (torque command) at n- 
th instant is expressed as follows: 

T e (n) — T e (n-1) + Kp COre(n) T Ki CO r e(n) (15) 


Where T e ( n ) is the torque output of the controller at the n- 
th instant, K p is the proportional gain constants and Ki is 
the integral gain constant. A limit of the torque command 
is imposed as 
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L e(n+l) - 


T T >T 

emax 5 e(n+l) emax 

"^emax ^~^e(n+l) e 


(16) 

The gains of PI controller shown here are selected by trial 
and error method, 


5. Results and Discussion 

The dynamic performance of the PMSM drive 
is tested for three different cases as mentioned 
below: 

i. Balanced Three Phase Supply 

ii. FOC 

iii. MRAS 

i. Balanced Three Phase Supply 

At first, the motor starts on no-load, after reaching to its 
rated speed of 1500rpm, a load torque of Tl= 2Nm is 
applied at t=0.1sec. 

In the first case, the dynamic performance of PMSM Drive 
under balanced three phase supply is modelled in 
MATLAB/ Simulink/ Sim Power System environment. 
Response of Stator currents, speed, torque, direct current, 
quadrature current and flux waveforms with three phase 
balance supply are plotted in Figure 5. In this case, there 
are lot of disturbances found in reaching steady state for 
speed and torque responses during starting. Whereas 
torque response has disturbance or lot of fluctuations when 
load is applied at t=0.1 sec. Response of torque and flux 
ripple waveforms with three phase balance supply are 
given in Figure 6. Since the input to the motor is a three 
phase balanced supply, it is observed that smooth response 
for torque and flux ripple waveforms as compared to the 
other two cases. 

ii. FOC 

In the second case, the dynamic performance of PMSM 
Drive with FOC is modeled. Response of Stator currents, 
speed, torque, direct current, quadrature current and flux 
waveforms for rated speed are plotted in Figure 7. 

In this case, there are few disturbances found in reaching 
steady state for speed and torque responses during starting 
as compared to the first case. Here the load torque Tl= 
2Nm at t=0.1sec; Ti = 1.5Nm at t=0.4sec and Ti = INm at 
t=0.5sec is applied and the dynamic performance of the 
drive is tested with rated constant speed 1500 rpm. 
Responses of torque and flux ripple waveforms with FOC 
for rated speed are given in Figure 8. 

Also the drive performance is tested considering constant 
rated load torque of Tl= 2Nm is applied at t=0sec with 
varying speed from rated speed; 1500 rpm at t=0 sec; 
lOOOrpm at t=0.1sec and 500 rpm at t=0.4sec. Response of 
Stator currents, speed, torque, direct current, quadrature 
current and flux waveforms for rated load torque are 
plotted in Figure 9. Responses of torque and flux ripple 
waveforms with FOC for rated load torque are given in 
Figure 10. The torque and flux ripples are quite high as 
compared to the first case. 
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iii. MRAS 

The same analysis is repeated for third case, the dynamic 
performance of PMSM Drive with MRAS is modelled. 
Response of Stator currents, speed, torque, direct current, 
quadrature current and flux waveforms for rated speed 
with MRAS are plotted in Figure 11. Responses of torque 
and flux ripple waveforms with MRAS for rated speed are 
given in Figure 12. Response of Stator currents, speed, 
torque, direct current, quadrature current and flux 
waveforms for rated load torque with MRAS are plotted in 
Figure 13. Responses of torque and flux ripple waveforms 
with FOC for rated load torque with MRAS are given in 
Figure 14. 

It is found that the drive performance is better in terms of 
speed and torque wave forms and also the torque and flux 
ripples are quite low as compared to the FOC. The 
parameters used for PMSM is given in Table 1. 
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Figure 5. Response of Stator currents, speed, torque, direct 
current, quadrature current andflux waveforms with three 
phase balance supply 



Figure 6. Response of torque and flux ripple waveforms 
with three phase balance supply 
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Figure 7. Response of Stator currents, speed, torque, 
direct current, quadrature current andflux waveforms 
with FOC for rated speed 1500 rpm 
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Figure 8. Responses of torque and flux ripple 
waveforms with FOC for rated speed . 
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Figure 9. Response of Stator currents, speed, torque, 
direct current, quadrature current and flux waveforms 
with FOC for rated load torque 2 Nm 




Figure lO.Responses of torque and flux ripple waveforms 
with FOC for rated load torque 
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Figure ll.Response of Stator currents, speed, torque, 
direct current, quadrature current andflux waveforms 
with MRAS for rated speed 1500 rpm 
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Figure 12. Figure 7. Responses of torque and flux ripple 
waveforms with FOC for rated speed . 
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Figure 13. Response of Stator currents, speed, torque, 
direct current, quadrature current and flux waveforms 
with MRAS for rated load torque 2 Nm 



Figure 14.Responses of torque and flux ripple waveforms 
with MRAS for rated load torque 


Table 1: Parameters of PMSM 


Parameter D e scription 

Value 

DC Voltage (v d O 

230 V 

Direct Axis Inductance^ L j 5 ) 

0.075 H 

Inductance 

0.095 H 

Rot or Flux C onst ant { Af) 

0-1540M 

Stator Resistance^ R 5 ) 

ISfl 

Numb er o f Poles [ p) 

4 

Moment of Inertia { J) 

C'.Sl 1 1 e SK *el : 

Viscous Fiiction Co-efficient (B) 

0.00 1 1 Nm 3 

Speed (N) 

1500rpm 


6. Conclusion 

Sensor less Control of PMSM with FOC using MRAS for 
speed estimation is carried out using MATLAB/ 
Simulink/Sim Power Systems is carried out in this paper. 
The output of the adaptation mechanism is the estimated 
speed quantity, which is used for the tuning in adjustable 
model and also for feedback. The stability of such closed 
loop estimator is achieved through Popov’s Hyper stability 
criterion. The method is simple and requires less 
computation. The drive is tested for three cases and also 
for dynamic conditions. It is found that with MRAS with 
PI Controller, the drive performance is smooth and torque 
and flux ripples are less. Further this work can be extended 
with Fuzzy Logic Controller along Multilevel Inverter 
topology.lt is further suggested that the proposed work 
further can be implemented with MRAS using fuzzy logic 
and FPGA Hardware implementation. 
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Abstract 

An optimal control problem with a linear state equation and vector control with the minimization of a 
quadratic criterion is considered. Furthermore, the terminal state is fixed with fixed horizon. The relaxation 
method is proposed coupled with the shooting technique to solve this problem. We rely on shooting method 
to find initial condition of adjoint vector p(0). The convergence of the numerical procedure is presented 
where numerical experiments are also presented. We solve the problem analytically then we compare the 
results. 

Keywords: Optimal control, Shooting method, Relaxation method. 


1 Introduction 

The theory of optimal control has been well developed for over forty years. With the advances of computer 
technique, optimal control is now widely used in multi-disciplinary applications such as biological systems, 
communication networks and socio-economic systems etc. As a result, more and more people will benefit 
greatly by learning to solve the optimal control problems numerically. Realizing such growing needs, books 
on optimal control put more weight on numerical methods. In retrospect, [11] was the first and the "classic" 
book for studying the theory as well as many interesting cases (time-optimal, fuel-optimal and linear quadratic 
regulator (LQR) problems). Necessary conditions for various systems were derived and explicit solutions were 
given when possible. 


In this work, we are interested in linear quadratic optimal control problems with a fixed terminal time. 
Moreover the state is subject to some constraints and a value of the final state is fixed and vector control. 
In order to use a numerical solution, we rewrite the equations of optimality from the Pontryagin’s minimum 
principle. Then we obtain a system in which, the differential equation describing the state is equipped with 
an initial and final conditions. Note also that, the costate equation derived from the Pontryagin’s principle is 
equipped with no initial condition or terminal one to be used algorithmically. In order to determine the initial 
condition of the costate we use in this study, the shooting method [5] coupled with the relaxation method (see 
[1] and [3]). Under suitable assumptions, we analyze the convergence of the considered iterative method. 


This study is organized as follows. Section 2 defines the problem and presents the shooting method. The Sec- 
tion 3 is devoted to the description of the relaxation method coupled with the shooting one. The convergence of 
the proposed method is presented in section 4, while Section 5 contains the results of the numerical experiments. 
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2 Statement of the Problem 

Let us consider the following problem: 

{ Found u G U a d such that , 

J(u) < J(v), V U G U ad , 

Where 

T 

J(U) = | J (( X - £ d ) 2 + ku 2 )dt ; 
o 

Xd corresponding at desired state, it is the control and k control the comparative weight given to the two 
components of the cost function, i.e. , the values of k are chosen in order to grant more weight to precision or 
energy expenditure to the final state xt fixed and minimizing the cost functional, and U a d is the closed convex 
domain. 

Under the following constraints: 


x\ = —bx i + ax 2 + ui, 
x 2 = ax i - bx 2 + u 2 , 
x(0) sa x 0 , 

k ^( T ) = 9, 


With G = (1, 2), 7 = 3, k, = (1, 0), g = 2 , T = 2 , a > 0, b > 0. 
Otherwise, We can reformulate the problem as follows: 


x\ = —bx i + ax 2 + i/i, 
x 2 = aaq - bx 2 +^ 2 , 

< x(0) = a; 0 , 

, ^i( T ) = 


The matrix form of the state system is written as follows: 

x\ \ ( b a \ ( x\ 




a —b 


X2 


+ 


1 0 
0 1 


ui 

U2 


with 


T 

1 1 

0 ) 


X-l - x ld \ 


( Ml 

\ 

/[(XI - X ld ,x 2 - X 2d ) 

0 

1 


x 2 - X 2d 

+ k(u 1 ,u 2 ) t 

u 2 


0 

V 

) 


\ J 

V 

J 


]dt. 


If we solve the previous problem by a numerical procedure, we have then to solve an algebraic-differential 
systemOs On one hand, the state equation describing the physical system is provided with an initial condition 
xq and the final condition Xf. On the other hand, the second equation corresponding to the costate one is 
formed with no initial condition and with no terminal condition. We will therefore use the shooting method to 
compute the value of p( 0 ). 


2.1 Shooting method 

The shooting method is used to obtain the value of p(0) necessary to the solution of the problem characterized 
by the Pontryaguin principle. If it is possible from the condition of minimization of the Hamiltonian to express 
the extremal control with respect to (x(t),p(t)), then the extremal system is a differential system of the form 
z(t) = F(t, z(t )), where z(t) = (x(t),p(t)). With a numerical integrator beginning from zq we obtain: z*° ~ z(U ), 
where the U are the discretized time generated by the numerical integrator. Note that in zo = (xo,Po), the value 
of xo is given by the initial condition of the problem and po denotes the value of p(0). Therefore, for different 
values of po, we will obtain the corresponding values of z *° . We are interested in the value of z 2 ^ ^ z{T), at the 
final time, where N is such that T = N.dt , dt being the time step, and z z ^ = In our case only x z £ 

are significant. Since the numerical results only depend on po they are denoted in the sequel x p ^ . We define 
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the mapping G as an implicit function which gives po by using a numerical integrator such that we obtain 
x ^ — xt = 0. Let 

G : M n -A M n such that G(po ) = x \ ^ — xt- 

Using G we define a nonlinear implicit system of n equations with n unknowns such that: 


G(po) = 0 . 


( 1 ) 


In order to solve (1), we will use the Newton’s method. The principle of the Newton’s method is described as 
follows: in the k th step, let be an approximation of the zero po of G; therefore po can be written po = Po + Ap§ 
and then: 


0 — G(po) — G(pq + A Po ) ~ G(Po) + 


dG 

dp 0 


{Po)-(PO - Po) + o(p 0 - Po), 


which leads to the solution of the following linear system: 


dG 

dpo 


( P k o)-(Po-p k 0 ) 


-G(p k 0 ), 


where §^(Po) is the Jacobian matrix of the application po — ^ G(po) computed when po = p§ ; note that the 
mapping po G(po) is not explicitly known but is known numerically. So we will use a method of numerical 
derivation based on the finite difference method. To avoid the computation of §p^(Po), if is sufficient to find 
an approximation of J^(p§). According to [4], we have two typical finite difference approximations which are 
frequently used: 


dG, 

dpoj 


( P k 0 ) 



or 

r)G 1 

JT^(Po) ~ T~[Gi(po + hije J ) - Gi(p 0 )}, 

OPoj tlij 

where the hij are the given discretization step parameters of the i th equation with respect to the j th variable 
and e k are the k th vector of the canonical basis; note that, classically, we can always choose the values of h^ 
equal to each other at a common value h. Let A ,j(p, h) be a finite difference approximation, then we have: 


Let 

J(p, h) = (A,j(p 0 , h)), 

which is an approximation of the Jacobian matrix; then the approximate Newton’s method can be written as 
follows: 


Po +1 =P k o- J(p k o, h k )-\G{p k 0 ), k = 0, 1, , 

The problem of convergence of this iterative process is ensured by using a result of the book of Ortega and 
Rheinbold [4]; indeed, if the discretization step h^ are small enough and tend to zero then the convergence is 
ensured. 


3 Numerical Algorithm 

To solve the problem, with and without constraints, we perform the coupling of the relaxation method (see [1] 
and [3]) with the shooting method [5], the latter method being intended to calculate p(0) necessary to solve 
the differential-algebraic system obtained by applying the Pontryagin minimum principle. The steps in this 
numerical method are summarized below: 

In the unconstrained case note that x = x ^ . 

Remark 1 Steps (2), (2) and (S)correspond to the relaxation method while the following steps correspond to 
the implementation of the shooting method. 

Definition 1 An invertible matrix A is an M-matrix if A -1 > 0 and a ij < 0 for i 7 ^ j. 

Remark 2 M-matrices have many main properties, in particular the spectral radius of the associated Jacobi 
matrix J = I — D~ x .A, where D is the diagonal of A, is lower than one; this property will be used in the sequel. 
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Algorithm 1 Numerical Algorithm 

<0> Initial Approximation 

Control u° for t G [0 , T], and of the costate variable p°(0), 
o — r i — 0 (where r is the number of iteration), 

<0 If | — u |> e (%where e defines the convergence threshold) do: 

- Determine the state variable x^ r \ by integration of the state equation: 


r <H = Ax + Bu ( < r \ 0 < £ < T, 

\ x(0) = x 0 , 

andx = Proj(x ), 

where Proj{.) is the projected operator, 

-and find the costate variable p ^ by solving: 


J dt 
\ P (r H 0), 

where p ^ (0) is computed by the shooting 
method, 

— Determine the control u^ r+1 ^\ 


(2) 


, 0 + 1 ) 


<— ( U ^ 




(3) 


o— Else | u^ r+1 ^ — \< e , 

- Determine the shooting function: 


G(p) = x (r) (T) - x f , 


- Solve the shooting equation by the Newton 
method and find the new value of p(0): 


p( r+1 )(0) <— p^ r \ 0) + correction , 


- r r + 1. 
§End 
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Proposition 1 If the following conditions are satisfied: 

- A is an M-matrix 

- k > ko > 0 
-p 2 (0)-p 2 (T) >0, 

then the algorithm for computing the numerical optimal control law, by the relaxation method coupled with the 
shooting method, converges whatever be the initial value of u° . 

4 Numerical Experiments 

The Hamiltonian of the system is given by: 

H(x(t),p(t),u(t),t) = i[(xi - Xidf + (x 2 - x 2 d) 2 + k(u\ +u%)\ 

+ (pi (t), P 2 (t)) ( ~ b H ) + (pi(t)Mt)) ( H ) 

= ^[(*1 - x l d) 2 + ( x 2 - X 2 d) 2 + k{u\+u\)\ +Pl(t)(-bx 1 + ax 2 + Ml) 

+ p 2 (t)((axi - bx 2 + u 2 ) +Pi(t)ui + p 2 (t)u 2 . 

On the state the optimality equations can be written as follows: 

f f ! & = % = -bxi+ax 2 +u 1 , 

< % = %F a =°*i-bx2 + u 2 , 

{ x(0) = x 0 , 

< ( = ~h>i+ a P2+xi-x ld , pi(T)=ai 

{ -^dt = = a Pi - b P2 + x 2 - x 2d, P 2 {T) = a 2 

PM- = 0 = kui + P i(t). 

M = 0 = +P 2 X). 


where a^,i = 1, 2 is found below. Adjunct on the previous system transversality condition: 
3 k\ A 0 such that: 


¥>(0, x 0 , T, x(T), ko, k\) = (^i(T,x(T)|A;i), 


with 


4>i (T,x(T)) = kx(T) - g , 


<fi(0,xo,T,x(T), ki) = (kx(T) -g\ki). 


P(T) 


dcp 

dx(T) 


(0, xq, T,x(T), k\) = -K l ki, 


( Pl (T), P2 (T)) = -(h0)h = 


After, the values of oti,i = 1,2 are defined, Hamilton-Pontriaguine equation’s are written as follows: 


{ x\ = —bx 1 + ax 2 + u\, 
x 2 = ax 1 - bx 2 + u 2 , 
x(0) = x 0 , 

< f -pi = -bpi + ap 2 + xi - x ld , pi(T) = -ki 
\ ~P 2 = api - bp 2 + x 2 - x 2d , p 2 (T) = 0 
0 = kui + P i(t) 

0 = ku 2 +p 2 (t) 
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Then, the matrix form is given as follows: 


( X 1 ^ 


/ ~b 

a 

0 

0 

0 

0 

0 



f X! 

\ 

X2 


a 

-b 

0 

0 

1 

0 

0 



X2 


-Pi 


-1 

0 

-b 

a 

0 

0 

0 



Pi 


~P2 

= 

0 

-1 

a 

-b 

0 

0 

0 



P2 


0 


0 

0 

1 

0 

k 

0 

0 



U\ 


0 


0 

0 

0 

1 

0 

k 

0 





V Pi( T ) ) 


V o 

0 

0 

0 

0 

0 

-1 

) 


\ ki 

) 


We solve iteratively the following problem of the fixed point: 


x\ = —bx i -f ax 2 + ui, 
x 2 = ax i - bx 2 + u 2 , 

^i(O) = 0,aj 2 (0) = 0, 
pi = bpi - ap 2 -xi+ xi d , 

P 2 = -api + bp 2 -x 2 + X 2 d, 

U 2 = ~^, 


Pi(T) = -k 1 
P2(T)= 0, 


5 Analytical Solution 

Let us derivation method at the level equation. We obtain: 

Pi = bp i - ap 2 -xi+ x ld , 


Pi = bpi — ap2 — x \ , 

Pi = b(bpi - ap 2 - xi + xi d ) - a(-api + bp 2 - x 2 + x 2d ) ~ {-bx i + ax 2 - y ), 
Pi = b 2 pi - abp 2 - bx i + bxi d + a 2 pi - abp 2 + ax 2 - ax 2d + bx i - ax 2 + 


Then, 


In the same way: 


pi = (a 2 + b 2 + ^)pi-2abp 2 + bxi d - ax 2d , 
pi = (a 2 + b 2 + t)p! - 2o6p 2 + bxid - ax 2d - 


P 2 = ~api + bp 2 -x 2 + x 2d 


P 2 = -api+bp 2 -x 2 , 


p 2 

p 2 = -a{bpx - ap 2 - xi + xid) + b{-ap 1 +bp 2 - x 2 + x 2d ) ~ {axi -bx 2 - 


p 2 = -abpi + a 2 p 2 + ax i - ax u - abpi + b 2 p 2 - bx 2 + bx 2d - ax i + bx 2 + 


p 2 = -2abpi + ( a 2 + b 2 + -)p 2 - ax u + bx 2d , 
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Then 


i >2 = (a + b 2 + -)p 2 - 2abpi - ax ld + bx 2 d 


derive double the equation (4), we obtain: 


(4) 


( 4 ) 

Pi 


(o + b + -)pi - 2abp 2 , 


p = (a 2 + b 2 + -)pi - 2ab(-2abpi + (a 2 + b 2 + -)p 2 - axid + bx 2 d), 


Then, 


(4) leads: 


p i 4) = (a 2 + b 2 ^)pi + 4a 2 b 2 pi - 2ab(a 2 + b 2 + ^)p 2 + 2a 2 bxi d - 2ab 2 x 2 d), 


p J 4) = (a 2 + b 2 + ^)pi + &a 2 b 2 p\ - 2 ab(a 2 + b 2 + t) + 2a 2 bxi d - 2ab 2 x 2d , 


2abp 2 = ( a 2 + b 2 )pi - pi + bx u - ax 2d , 


By injecting (5) in (5), we obtain: 


(5) 


Pi 1) = (a 2 + b 2 + y)pi + 4a 2 b 2 pi - ( a 2 + b 2 + y)((a 2 + b 2 + y)pi - pi + bx u - ax 2d ) + 2a 2 bxi d - 2ab 2 x 2d , 


p j 4 ' 1 = (a 2 + b 2 + \)pi + A.a 2 b 2 pi - ((a 2 + b 2 ) 2 + y(a 2 + b 2 ))pi + (a 2 + b 2 + ypi 

k k k 

- b(a 2 +b 2 + y)xi d + a(a 2 + b 2 + \)x 2d + 2a 2 bxi d - 2 ab 2 x 2d , 

k k 


p ^ = (2a 2 + 2b 2 + |)pi - ((a 2 + b 2 + ^) 2 - 4a 2 6 2 )pi + (-a 2 6 - b 3 — ^ + 2a 2 b)xi d 

a 


+ (u 3 + ab 2 + — — 2 ab 2 )x 2 di 


Then, 


p^ 4) - (2a 2 + 26 2 + |)pi - ((a 2 + 6 2 + t) 2 - 4a 2 6 2 )pi = (-6 3 - ^ + a 2 b)xi d + (a 3 - a6 2 + ^)x 2 d , 

The characteristic equation corresponding in equation (6) written as: 

C 4 - 2 (a 2 + b 2 + l)C 2 - (-{a 2 + b 2 + \) 2 + 4 a 2 b 2 ) = 0 , 
k k 

A = 4 a 2 b 2 

The root of the characteristic equation are given by: 

C\ = (a 2 + b 2 + y) — 2 ab, 


(6) 


= ( a 2 + b 2 + — ) + 2 ab. 
k 


Then 


Pi{t) = A e Clt + /3e Clt + fie C2t + ae ° 2t + v. 


(7) 
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Determine z^, we have: 

pi(i) = \C ie Clt - + /iC 2 e C2t - aC 2 e- G2t , 

i>i(i) = A Cfe Clt + pC^e~ Glt + pCle° 2t + aCle~ C2 \ 

p (3 \t) = A Cfe Clt - /3Cle~ Clt + /iC%e C2t - aC%e~ C2t , 

p[ A \t) = XCfe Clt + pC\e~ Clt + pC\e C2t + aC\e~ C2t . 

We substitute on (6), and we obtain: 

A [Cf - {-(a 2 + b 2 + t) 2 + 4 a 2 b 2 ) - 2{a 2 + b 2 + i)C?]e Clt + /?[C a 4 - (-( a 2 + 6 2 + t) 2 + 4aV) 

- 2 (a 2 + b 2 + t)C 2 ]e- Clt + M [C 2 4 - (-(a 2 + & 2 + t) 2 + 4 aV) - 2(a 2 + b 2 + t )C|]e C2t + a[C 4 

- (-(a 2 + 6 2 + b 2 + 4a 2 b 2 ) - 2(a 2 + 6 2 + bcfje -02 * + (-(a 2 + 6 2 + t) 2 + 4 a 2 b 2 )u 

rv rv rv 

= (~b 3 - y + a 2 b)x ld + (a 3 - ab 2 + y)x 2 d, 
k k 

By identification, we obtain: 

_ (-& 3 - | + a?b)x i d + (a 3 - ab 2 + %)x 2d 
4 a 2 b 2 — ( a 2 + b 2 + j:) 2 

We deduce P 2 (t) of (5), we obtain: 

P 2 (t) = T[( a 2 + b 2 + k)pi -pi + bx ld - ax 2d 2ab], 


p 2 (t) = — ——r— ^ [Ae Clt + /3e~ Clt + pe° 2t + ae~ C2t + v] - — [\C 2 e G ^ + pC 2 e~ c ^ + pC 2 e G2t 

ZCLO ZCLO 

+ aC 2 e~ C2t } + -^xi d - ^x 2 d , 


P2(t) = A[- 


i 2 + b 2 + l- C 2 


+ p[ 

+ ( 


2 ab 

a 2 + b 2 + i — C 2 


* -]e Clt + /?[- 


b 2 + \~C 2 _ 


2 ab 
a 2 +b 2 + 1 


- — — ]e C2t + a[- 


2 ab 
b 2 + h-C 2 


Cit 


2ab 


l] e -c 2 t 


2ar k)p+ h xu - : k x ^ 


Then, 


p 2 (t) = \e Clt + fie Clt -pe° 2t -ae ° 2t + ( ° + 2 ^ + fc )^ + Tx id - ^x 2 d- (9) 

Also, given the following equation: 

a-'i = &Pi - ap 2 - pi +xw, 
x 2 = -api + bp 2 -p 2 + x 2d . 

We obtain the following results: 


Xi(t) 


+ 


6[Ae Cl ‘ + /3e~ Clt + pe° 2t + ae~ C2t + v] - a[Ae Clt + /3e~ Clt + pe° 2t + ae~° 2t 
+ 2afe + ~ V + ‘L' Xld ~ h) X2d ^ ~ ^ Cie Clt - /3Cie~ Clt + pC 2 e° 2t - aC 2 e~ C2t ] 


+ Xi d 
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x\ (t) = A(b — a — C\)e Clt + /3(b — a + C\)e Clt + fi(b + a — C 2 )e° 2t + ot(b + a + C 2 )e ° 2t 

r, (& 2 + b 2 + ^) 1 a 

+ [b ~ b \u - 2 .r„, I -x 2d I 


Then, 


xi(t) = A(6 — a — Ci)e Clt + /3(6 — a + Ci)e Clt + fi{b + a — C 2 )e° 2t + a(6 + a + C^e C2t + 


b 2 - a 2 -h 


2 b 


1 a 

+ 2* ld + 26*“- 


( 10 ) 


x 2 (7) = -a[\e Clt + f3e~ Clt + /je C2t + ae~ C2t + u\ + b[\e Clt 

+ (b 2 - a 2 - t)2% + ^x ld + 


Then, 


+ /j,e C2t + ae~° 2t 

7^ x 2 d ~ [\Cie Glt 


+ /3e~ Clt + ne C2t + ae~ C2t + 
/?C 2 e _Clt + n C 2 e G2t - aC 2 e~ G2t ] + x 2 d, 


X2 (t) 


= A (b — a 
b 2 -a 

+ 


- Ci)e Glt +/3(b-a + Ci)e- G ^ + n(-b - 

•• r !. b 1 

^ + -x la + -x„. 


■ a 


- C 2 )e G2t + a(-a -b- C 2 )e~ G2t 


+ 2a P+ Ya XU+ 2 Xm - 

The constants are determined by the following li: 

•n(0) = o,x 2 (0) = o, 

xi(T) = 2, 

m (T) = -jfei.cm = 0. 


( 11 ) 


limits conditions: 


■i V 1 ) = 
h(T) = -k 1 ,p 2 (T) = 0. 

Ve solve the linear system numerically: 

0 


= X(b- 


a + Ci) + fi(b + a — C 2 ) + cx(b + a + C 2 ) H 25 — ~ 

0 = A(6 
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: to following system: 
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Ci)c^ T 
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V 2a6 2a 
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Let matrix form: 
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e c lT 
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e~ ClT 
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(— b - a + C2)e G2T 

e c 2 r 

_ e c 2 T 


b + a + C 2 


-b-a-C 2 

(— b — a — C 2 )e~° 2T 

e -c 2 T 


-e-W 


0 \ 
0 
0 
1 

0 / 


A ^ 


p 



= 

a 


\ki ) 



fe 2 -° 2 -j 


„2_ 1 
fc 

2 b 

(b_ ^C±t 
[ 

0 b —a — 


2b 


-( 


a 2 +fe 2 + l 
2 ab 


-V - \xid - ^x 2d ^ 

)y-^xid-\x 2d 

-V - 2 X\ d - j^x 2d 
—v 

)v-^x ld + ±x 2d / 


5.0.1 Comparison of the Two Approaches 

The numerical algorithm is implemented in Matlab. Particularly, we used the function ode 45 et f solve. In 
this example, the method converge independently with initial point ^(0) = (0.0536,0.1072). Numerical 
experience is determined by values a = 1 , b = 3 ct k = 2. Let us deduce that analytical solution and numerical 
solution correspond perfectly. For different values of &, The performance of the numerical procedure is 


( a ) 


(b) 




Figure 1: (a) u\ numerical solution (b) v\ true solution 
(a) (b) 




Figure 2: (a) u 2 numerical solution (b) v 2 true solution 


summarized in the following table. 


k 

time C.P.U 

Iteration Number 

0.5 

1.3089 

12 

1 

2.3619 

8 

1.5 

2.7016 

7 

2 

3.7995 

7 

2.5 

5.7479 

6 


Note that the algorithm converges quickly, to a very small number of iterations required to reach convergence. 
In addition, the computing time used is very low. 
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(a) 


(b) 



Figure 3 : (a) x\ numerical solution (b) y\ true solution 
(a) (b) 
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Figure 4 : (a) x 2 numerical solution (b) 7/2 true solution 
(a) (b) 



Figure 5 : (a) p\ numerical solution (b) q\ true solution 
(a) (b) 



Figure 6 : (a) P2 numerical solution (b) true solution 


6 Conclusions 

The goal of this study concerns the linear quadratic optimal control problems with a fixed terminal time. We 
have used the relaxation method coupled with the shooting method for solving the optimal control problem. 
The convergence of the procedure is ensured. The rate of convergence is high and the computation of time is 
fast. We concluded that relaxation method give a good results in quick time. 
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Abstract 


1 Introduction 


Converters are employed in wide range of applications 
to save energy and attain desirable voltage. Matrix 
converter can convert three phase AC input to three 
phase AC output with variable voltage amplitude and 
frequency directly. It can be used as a bidirectional 
power flow converter without any intermediate stor- 
age element. The objective of this proposed research 
is to minimize the harmonic losses to get maximum 
output voltage ratio, sinusoidal current, desired vari- 
able voltage amplitude and desired variable frequency. 
Pulse width modulation and space vector modulation 
algorithm control the input and output voltage and 
frequency independently. In this research work, 
pulse width modulation and space vector modulation 
based matrix converters will be designed to attain 
voltage ratio up to 1 and to reduce switching losses 
so that total harmonic distortion could be minimized 
with sinusoidal waveforms of desired amplitude and 
frequency. Comparisons will be drawn on the basis 
of harmonic contents present in the desired output 
waveform by changing the characteristics of input 
waveform. Simulation environment will be created in 
Mat lab. 1 

Keywords: Pulse Width Modulation, Space Vec- 
tor Modulation, Matrix Converters, Matlab 

N omenclat ure 


PWM 

Pulse Width Modulation 

SVM 

Space Vector Modulation 

AC 

Alternating Current 

DC 

Direct Current 

BJT 

Bipolar Junction Transistor 

MOSFET 

sistor 

Metal Oxide Silicon Field Effect Tran- 

IGBT 

sistor 

Insulated Gate Bipolar Junction Tran- 


Mhis study has been implemented at Department of Elec- 
trical Engineering, University of Lahore, Pakistan 


The basic component of matrix converter is bidi- 
rectional switch which are forcibly commutated. 
Switches are used in controlled manner to generate 
a variable frequency output. There is no dc link be- 
tween two different voltage sources and also there is 
no need of using energy storing active elements link 
capacitor. The switches are used in controlled way 
due to which the operation at high frequency is pos- 
sible. There are many type of matrix converters in 
which the components used can be reduced according 
to the requirement and thus the modulation technique 
is also less complex [22, 9, 20, 18, 15]. 

Pulses of variable widths are generated using pulse 
width modulation that represent the amplitude of an 
analogue signal. The main use of pulse width mod- 
ulation is that the power supplied to the electrical 
loads especially of rotating loads like motors which 
have some inertia, can be controlled. Thus, motors 
can be made to operate at variable speed by using 
pulse width modulation. 

To control the pulse width modulation, another 
algorithm space vector pulse width modulation is em- 
ployed [19]. It is used to drive a motor at variable 
speed by generating three phase alternating currents 
waveform. There are many methods and variations 
of space vector pulse width modulation, each having 
different quality and different requirement of compu- 
tation. Due to rapid switching actions, the harmonics 
are generated. The elimination of these harmonics is 
one of the basic purposes of using space vector pulse 
width modulation [30]. 

Non linear loads generate current harmonics and 
voltage harmonics in a system. High order frequen- 
cies are major cause of problems in quality of power 
supplied. Harmonics result in increased heating of the 
equipments. It is highly desirable to reduce the har- 
monics in any power system. Low order harmonics 
which are undesirable, can be eliminated by select- 
ing most preferable harmonics elimination technique 
[21, 24]. The operation of electrical system is highly 
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effected due to generation of harmonics. A lot of at- 
tention in any system is paid to the reduction of har- 
monics [11, 27]. The operations of a systems have 
standards to maintain and to limit up to which har- 
monics can be added to the system [1, 2, 3]. By using 
multiple winding transformer in multiple connection 
at input side of a system, current harmonics can be re- 
duced. Harmonics can also be removed using different 
active and passive filter topologies [5, 28, 4]. 

The matrix converter topologies are more attrac- 
tive after the development of power switches like bipo- 
lar junction transistors (BJTs), metal oxide silicon 
field effect transistors (MOSFETs) and insulated gate 
bipolar junction transistors (IGBTs) [12]. The first 
time development of matrix converter was in 1980 
[31, 32]. Venturini and Alesia in 1980, presented the 
converter using bidirectional switches and they in- 
troduced the name of matrix converter for the first 
time. The modulation technique introduced by them 
is known as direct transfer function approach. In or- 
der to obtain the desired output voltages, the input 
voltages are multiplied by the modulation matrix [23] 
[ 8 ], 

After them in 1983, Rodriguez introduced a differ- 
ent technique for matrix converter operation introduc- 
ing fictitious DC link [29]. According to his technique, 
the output was switched between the most positive in- 
put line and the most negative input line. This tech- 
nique is known as indirect transfer function approach 
[26]. 

The method for the control of Matrix Converters 
was introduced by Braun in 1983 [6] and by Kast- 
ner and Rodriguez in 1985 [14]. The first paper to 
provide the solution of modulation scheme in matrix 
converters using space modulation was published in 
1989 by Huber [10]. Schauder and Neft proved this 
fact experimentally in 1992 that high quality of input 
current and output current can be obtained by using 
matrix converter which comprises of only bidirectional 
switches [25]. The switching should be performed in 
such a manner that there should be no high spikes of 
voltages and current which can damage the switches 
made up of semi-conductor materials. This fact af- 
fected the interests of using matrix converters. This 
issue has been solved due to the development of micro- 
processors based controllers and multistage switching 
and commutation techniques. 

Implementation and operation of bidirectional 
switches [17, 7] under normal condition with active 
filters [16, 33] are main themes of researches. Lanka 
Parampil, Ushakomari and Nisha removed the har- 
monics in three phase inverters using space vector 
pulse width modulation [13]. 

In this research work, three phase input and three 
phase output will be compared and modulated accord- 
ing to requirement by a modulation block. The mod- 
ulated voltages will be fed to a matrix converter which 
comprises IGBTs whose operation will be controlled 


by applying a PWM and SVM at the gates of the IG- 
BTs. The output from the matrix converter will be fed 
to different loads like synchronous motor, induction 
motor, resistive load and RL load whose harmonics 
will be removed. 

The remainder of the paper is organized as follows: 
Section II focuses on description and mathematical 
framework of matrix converter, pulse width modu- 
lation and space vector modulation Section III em- 
phasizes on mathematical formulation of modulation 
techniques. Section IV covers the simulation results 
and Matlab/Simulink models. This research work and 
its simulation results are concluded in section V. 

2 Mathematical Formulation of 
Research Problem 

2.1 Matrix Converter 

A device that changes the frequency of input supply 
is called matrix converter. The basic components of 
matrix converter is bidirectional switches which allows 
any output phase to be connected to any of the three 
input phase for a given length of line and any desired 
frequency can be generated. The main advantages of 
the matrix converter are as follow: 

• It eliminates the need of large reactive compo- 
nents which are used in AC-DC-AC type invert- 
ers. 

• The use of bidirectional switches allows the en- 
ergy to be regenerated back to the supply or 
grid. 

• This allows the input current waveforms to be si- 
nusoidal and suitable switching can be provided 
the unity input power factor. 

• Due to absence of large reactive elements like 
capacitors and inductors, the size is small. 

• Matrix converter provides more power to weight 
ratio. 

• DC link in back to back inverters is absent. 

• The number of switches used are also less. 

2.2 Different Types of Matrix Con- 
verters 

2.2.1 AC-DC-AC Converters 

These types of converters are known as dc-link con- 
verters. There are two types of dc-link converters, 
which are voltage source inverters and current source 
inverters. In voltage source inverters, the rectifier 
stage is realized using diodes bridge while the dc link 
circuit consists of a shunt capacitor. In current source 
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inverter the rectifier stage is realized using phase con- 
trolled switching bridge and he dc link circuit consists 
of serried inductor. An AC-DC-AC converter with 
bidirectional power flow can be realized using pwm 
rectifier and then a pwm inverter to the dc link. En- 
ergy storing element that is common to both sides 
forms the dc-link between them. 

2.2.2 Cyclo-Converters 

There is no dc link in cyclo-converters. The cyclo- 
converter generates the output of variable frequency. 
The output generated is almost sinusoidal because 
the segments of the input waveform are transferred 
to the output side by switching. SCRaAZs are usu- 
ally used in cyclo-converters. The output frequency 
of the cyclo- converter can never be greater than the 
input frequency. 

2.2.3 Matrix Converter 

Matrix Converter converts directly AC to AC with- 
out using any dc link. This increases reliability and 
stability of the system. Bidirectional switches using 
IGBTaAZs are used commonly in matrix converters. 
Depending on the number and types of components 
used, the matrix converters are further classified into 
many types. 

a. Sparse Matrix Converter 

In this type of matrix converter, the numbers of re- 
quired switches are less, so the complexity of the gate 
drive circuit is reduced. The function is identical to 
direct matrix converter. 18 diodes and 15 switches are 
required for the sparse type matrix converter. 

b. Very Sparse Matrix Converter 

In this type, the number and diodes are increased and 
correspondingly the number of switches is reduced as 
compared to the sparse matrix converter. Though 
gate drive complexity is reduced but due to increase 
in the number of diodes, the conduction losses are in- 
creased. 

c. Ultra Sparse Matrix Converter 

Only unidirectional switches are used in the input 
stage of ultra sparse matrix converter. So, these types 
of matrix converters are used for the variable speed 
drives which are of low dynamics. The topology of 
ultra sparse matrix converter introduces phase dis- 
placement between input voltages and input currents. 
Only 12 diodes and 9 switching devices are required 
for this type of matrix converter. 

d. Hybrid Matrix Converters 

The hybrid matrix converter converts AC/DC/AC 
but does not use any dc link or reactive elements like 
capacitor or inductor. The hybrid matrix converters 
are further classified into two types depending on their 
operation. If hybrid matrix converter converts both 
voltage and current commutation in same stage then 
this is called hybrid direct matrix converters. If cur- 
rent and voltages are converted in different steps then 


these are called hybrid indirect matrix converter. 


2.3 Basics of Matrix Converter: 


The matrix converter is direct AC- AC converter for 
converting one frequency AC supply to another fre- 
quency AC supply without involving an intermediate 
DC link capacitor. 

It has a three phase input supply. The three phase 
output voltages obtained are V a , Vb and V c . There 
are nine bidirectional switches from Sn to S 33 which 
represent a nine matrix components mathematically. 

The three phase matrix converter converts the 
three phase input of given amplitude (Vi) and fre- 
quency (fi) to three phase output of a fixed amplitude 
(y o ) and frequency (/ 0 ). Any desirable output fre- 
quency can be achieved by this converter. The three 
phase input voltages of the converter and required out 
put voltages are given by (1) and (2). 


'Va' 


COs(c Oit) 

V b 

= Vi 

COs(c dit + 

Vc_ 


COs(c Oit + 

'Va 


0 

0 

w 

£ 

v b 

= V 0 

COS (iJit + 

V c 


COS (iJit + 


The input and output voltage are related to each 
other according to the following matrix equations as 
shown in (3). Here M\j is the duty cycle of switch Sn 
and so on. 


Va(ty 


Mn Mi 2 Mi 3 


~V A (t) 

Vb(t) 

= 

M21 M2 2 M2 3 


v B (t) 

Vc(t)\ 


M31 M3 2 M3 3^ 


[Vc(t)\ 


Normally, the Matrix Converter is fed by a volt- 
age source and switching should be performed in such 
a way that none of the input terminals are short cir- 
cuited. The matrix converter applications are mostly 
used to drive inductive loads, so the the switching 
should also be performed such that the output phase 
is never open circuit because it will cause discontinu- 
ous supply to load and also the inductive kicks will be 
caused. The load and source voltages and input and 
output currents can be expressed as vectors defined 
by (4) and (5). 


Vo = 

Va{t) 

v b {t) 

Vc(t)_ 

;V i = 

'Va (*)' 
Vs(t) 

Vc(t ) 

(4) 

i 0 — 

ia{t ) 
ib(t) 

5 — 

u(*) 

i B (t) 

(5) 


2.4 Pulse Width Modulation 

Pulse width modulation is a technique in which a mes- 
sage or a signal is encoded in such a way that it takes 
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the form of a pulsating signal. This technique is used 
to encode any information that can be used for trans- 
mission. One of the main use of PWM is the con- 
trol of power that can be supplied to the load. By 
constantly turning on and off the switching device be- 
tween load and the source, the value of voltage re- 
quired is achieved. This phenomenon is carried out at 
high switching frequency. By varying the duty cycle 
of the PWM signal, the amount of power supplied to 
the load is varied. Due to constantly on and off the 
switching device the desired output waveform will not 
be smooth. So, in order to keep the output waveform 
smooth, the switching frequency should be as high as 
possible. 

The switching frequency of a PWM signal is very 
high which enables the power electronic switching de- 
vices to be saturated hardly. So, between on and off 
state of the switching the transition interval is very 
short and hence, switching losses are also less. Dur- 
ing off state of a controlled switch, the is no flow of 
current and during on state, the forward voltage drop 
is almost zero. So, using PWM signal for switching 
the switching losses are almost zero. 

Depending on the requirement, width of the pulse 
is modulated. The term duty cycle is defined as the 
ratio of on time of signal to the total time period of 
the signal as given by (6): 


given by (7). 




/ COs(u>it) \ 



V B (t) 

= 

cos(uJit + ?f) 


(7) 

\UV)J 


\cos(wit + 




The switching is performed after a sequence of 
time. This time T s which is reciprocal to the switching 
frequency f s . The switching time in term of switches 
shown is defined in (8) 

= tAa, T t Ba T tea 
= tAb + tsb + tcb 
= tAc + ^Bc + tCc 

- i 

The switching frequency f s is constant which 
means that the length of each sequence is same. The 
values that are to be received at the output side is say 
V a (t),Vb(t) and V c (t) which should be displaced from 
each other by 120 degrees. The three phase output 
voltages in terms of switching time are given by (9) 
(10) and (11) 


Duty Cycle (D) = t on /t on + t off (6) 

Duty cycle is represented in percentage like 50%, 
means on for half of the time and off for half of the 
time, 100% duty cycle defines fully on. 

The steady state operation of a dc to dc type con- 
verter is made possible when the reference signal used 
for PWM generation is constant and does not vary. 
In dc to dc converters, the reference is assumed to be 
a dc value. In AC to DC or DC to AC type convert- 
ers, the reference is assumed to contain the fundamen- 
tal frequency component of desired output frequency. 
Similarly in case of multiple or 3 phase converters, 
the sinusoidal reference signals are shifted by desired 
amount of phase shift as required in the output volt- 
ages. In three phase inverters, the reference signals are 
shifted by zero degree, 120 degrees and 240 degrees to 
generate the desired switching. In the generation of 
switching sequences, reference signals may also con- 
tain harmonics mostly 3rd harmonics which are self 
introduced so that utilization of dc voltages can be 
increased. 


2.4.1 Modulation algorithm for PWM based 
Matrix Converter 


Assuming that the switches are ideal, i.e., no losses 
occurs in switches. The three phase input supply is 


V a (t) = VimCOS^Cdt') ~~ + V im COs(ujt + 

+ cos(w£+y)^ (9) 

Vb(t) = Vi m cos{u)t) + V im cos(ut + 

+ cos(wt+y)^ (10) 

V c (t) = V im cos(u)t + V im cos(u)t + 

-L s 4 1 s 

+ COs(wi+y)^ (11) 

Now, if the input frequency is increased from uji to 
uj 0 then a modulating frequency uj m has to be added 
in the input frequency given by (12). 


U) Q — T ^ m (12) 

The output voltage V a (t) has zero degree shift in 
its phase and switches tBa and tea determine the 
V a (t) by connecting to first, second and third phase re- 
spectively. Thus to maintain zero degree shift in phase 
V a (t ) , the switching sequence tsa and tea should be 
retarded by 120 degrees and 240 degrees respectively. 
And same is for other two phases. The switching se- 
quences thus formed are given by (13) (14) (15) (16) 
(17) (18) (19) (20) and (21). 
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V a at 0° 


2.5 Space Vector Modulation 


tAa = y (1 + 2qcos(uj m t + 6)) (13) 

T 2,tt 

tBa = y (1 + 2qcos(co m t + 9 - —)) (14) 

T 4? r 

tea = y {1 + 2qcos(u m t + 6 - —)) (15) 

Vo at 120° 

T 47r 

tAb = y (1 + 2qcos(u m t + 9 - —)) (16) 

tBb = y (1 + 2 qcos(LO m t + 9)) (17) 

T 2ir 

tcb = y (1 + 2qco s (u m t + 9 - y)) (18) 

V 0 at 240° 

T 27T 

tAc = y (1 + 2qcos(u m t + 9 - —)) (19) 

T 47 r 

^Bc = y (1 + 2gcos(w m t + 6» - y)) (20) 

tec = y (1 + 2qcos(oj m t + 9)) (21) 


The generated three phase voltages vary from each 
other by 120 degrees with the frequency same as of 
reference signal. The reference signal can be varied 
by varying the time period T s as T s = l// s , where 
f s is switching frequency. The reference signal can 
be generated from a three phase using d — q or a — 
(3 — 7 transformation. Various combinations exist for 
selecting the switching sequence of switches but each 
strategy has its own switching losses as shown in figure 
1 . 
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Figure 1: This table shows combination of switching 
sequence. 


If these values of the time sequence of each switch 
is substituted in the equation of the desired output 
voltages as defined before, then the modulation matrix 
is given by ( 22 ). 


M(t) = | 


cos(uj m t) 
cos(u m t - ^f) 
cos(uj m t - ff) 


1 1 1 
1 1 1 
1 1 1 

cos(cj m t - ^f) 
cos(u m t) 
cos(co m t - ^f) 


I 

“r - 3 

cos(bj m t - ^f) 
cos(uj m t - Hp 
cos(uj m t 

( 22 ) 



Figure 2: Vector representation of SVM signal 


with uj m = (cjo — uji). The output and input voltages 
are related to each other with the help of modulation 
matrix as given in (23). 

V 0 (t)=M(t).Vi(t) (23) 

So the output voltages are given by (24). 


V(t) 


1 1 
1 1 
1 1 


1 

1 

1 


I 2 q 

3 


COs{(jJ m t) cos(u m t — ^ L ) 
COsibJ^t — ^) COsijjJmt) 

cos\ui m t - ?f) cos(L0 m t - ^f)) 


cos(w m t - ^f) 
cos(u} m t - ?f) 
cos(i0 m t) 


X 


V m 

Vm(oj m t - ff) 

y m (u m t - ^f) 


(24) 


Two zero vectors and six active vectors can be rep- 
resented in the form of a hexagon in such a way that 
the length of each vector represents the magnitude of 
voltages as showm in figure 2. The rotation of vectors 
represents the angular speed of the system. As pre- 
vious discussed for balanced loads, space vectors are 
represented in 2-Dimensional figure. To minimize the 
time taken during toggling of states and to maintain 
the effective frequency of the switches used at min- 
imum, the SVM topology prefers switching between 
two nearest active states. Using this approach, in ev- 
ery step only one leg is affected. Using proper switch- 
ing, the effective output voltage to input voltage ratio 
can be increased to more than 90.6%. 

The six active vectors can be represented by (25). 

Vk = ^e^-^with k = (1.2....6) (25) 
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SVM vectors are represented graphically by equally 
dividing six sectors of a hexagon. The switching of 
only one inverter leg provides minimum switching fre- 
quency and optimum harmonic performance. Using, 
SVM we can define any reference vector as a com- 
bination of two active vectors at any instant of time. 
Suppose that in figure 2, V re f shown is in k sector and 
the active vectors adjacent to V re f are V& and Vk+%* 
So, the reduced equation is given by (26). 


V re f = Vk.Tk + Vfc+i.Tjfe + i (26) 

Where T s is switching period. If this equation is split 
in real and imaginary parts, then it is given by (27) 
and (28). 
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(27) 
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3 3 
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sm cos 


T k 

Tk+i 

k is determined by taking the argument of the refer- 
ence vector given by (29). 


(28) 


(Jfe — 1) * pi 


< arg 


< 


k * pi 


(29) 


V re f makes a circular path of radius V re f at an angular 
velocity in the complex plane. So, larger the radius 
of the trajectory circle, largest is the magnitude of 
voltage that can be achieved. This circle is tangential 
to the mid points of the lines connecting the ends 
of the active state vectors. The maximum achievable 
phase voltages are given by (30). 


\Kef\ maic = \v d ^ = ±f d (30) 

Following the definition of modulation index intro- 
duced in above, the corresponding maximum modu- 
lation index is given by (31). 


TTlmax cont 


\V T 


ref |maa? 


73^ 


Vmax,sixstep _ r -Vd 2\^3 


= 0.906 
(31) 

With the definition of modulation index the compu- 
tation of the inverter switching times does not re- 
quire the knowledge of the adopted DC-link voltage 
but depends only on the desired modulation index. 


Only space vector topology provides the best suitable 
switching sequence according to the type of loads like 
pulsating load, non-linear load and static load etc. 
The advantages of SVM are: 

• THD of the output voltage is low 

• DC bus utilization in case of SVM is at least 
15% more than the PWM technique used 

• SVM offers low peak currents in controlled 
switches as compared to PWM. 

• Higher performance, efficiency and reliability 
is achieved using SVM as compared to PWM 
based inverters of similar type. 
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Figure 3: This table shows that 27 different combi- 
nations of switches are obtained using 9 bidirectional 
switches. Keeping in view the modulation constraints 
only 21 switching combinations are useful. The last 
three switching combinations provide zero vectors. 
The other six combinations are not useful because 
these combinations can not provide the reference vec- 
tors. 


3 Modulation Techniques 

The relations between input and output voltages are 
related to the states of the nine bidirectional switches 
with the condition that 0 < Mij <1 where = 
1,2,3. The variable represents the duty cycle of 
9 bidirectional switches [31, 32]. The duty cycle must 
satisfy the following equations (32) (33) and (34). 


Mu + Mi 2 + Mi 3 : 1 

M21 + M22 + M23 = 1 
M31 + M3 2 + M33 7 1 


(32) 

(33) 


(34) 
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Venturini proposed the method whose modulation 
function gives maximum value of voltage transfer ratio 
of 0.5 which is very low and is given by (35). 

1 27 r 

M i:j = - cos[6 0 - (j - 1) —] 

2 q .277 

+ — [cosa 0 - (i - 1)— J 

2t r 

* cos[6 e — (j — 1)— ] (35) 

This method was modified by optimum method to in- 
crease the voltage transfer ratio to 0.866 and its mod- 
ulation function is given by (36). 

1 27 r 

M ij = g{l + 2g[cos6 h -(j-l)y] 

* [cos(a 0 — (i — 1) y) 

- \ cos(3a 0 ) + (^\/3) cos(36/j)] 

6 2 

2 27 r 

- (~V3)qcos[(4:b h ~ (j ~ l)y) 

27 r 

- cos(2&h + (j — 1)— )]} (36) 


4 Simulation and Results 


The first solution is given by (37). 


M 1 (t) = \ 


1 

1 

1 


1 1 
1 1 
1 1 


cos(cu m t) cos(cu m t — ^) 
COs(cJ m £ — ) COs(cJ m £) 
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+ 2g 


COS^mt - 4p) 
COS^mt - ?f) 
COS^mt) 

(37) 


with cu m = (cuo — cJi). This yields pi = </>o, i.e. the in- 
put phase displacement is the same as the load phase 
displacement. The alternative solution is given by 
(38). 


M 2 (t) = | 


1 

1 

1 


1 1 
1 1 
1 1 


COs(cJ m t) C08(cd m t — 4^) 
C08(cd m t — 4p) C08(cd m t — ) 

C08(cd m t — 4^) C08(cd m t) 


+ 2q 


cos(uj rn t - 4p) 
cos(cd m t) 
cos(cu rn t - ?f)_ 

(38) 


with uj m = — (cuo + tdi). This yields pi = - 0 O , i.e. 
the input phase displacement is the reverse of the load 
phase displacement. Combining the two solutions pro- 
vides the means for input displacement factor control 
given by (39). 

[M(t)] = a 1 [M 1 (t)]+a 2 [M 2 (t)] (39) 


Matlab/ Simulink is used for simulation and results. 
A complete step by step process is shown in a figure 4. 
A three phase source is required whose amplitude and 
frequency can be varied. Modulation block also mod- 
ulates the switching intervals for the switches which 
in our simulation are insulated gate bipolar junction 
transistors (IGBTs). 




Figure 5: Gating Signal Generator 

Combined solution allows input displacement fac- 
tor control. This can be explained by taking the ex- 
ample of an inductive load. If: 

• al = a2 : input is resistive (unity displacement 
factor) 


Figure 4: Flow Chart for the Implementation of PWM 
and SVM based Matrix Converters 

The switching intervals are defined by using the 
solution as provided by the Venturini Method [31, 32]. 
Using Venturini method, two solutions are obtained. 


• al > a2 : input is inductive (lagging displace- 
ment factor) 

• al < a2 : input is capacitive (leading displace- 
ment factor) 

The output voltages are thus produced by us- 
ing these switching intervals and the input voltages. 
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Since, the transfer ratio is 0.5, so the output voltages 
cannot be increased by the half value of input voltages. 
Gating signal generator is shown in figure 5. The six 
signals of the switching intervals are compared with 
the repeating sequence taken as reference to produce 
the gate signals for the IGBTs. This is pulse width 
modulation. Matrix converter and space vector mod- 
ulation based matrix converters are shown in figures 
6 and 7 respectively. 




Figure 8: Circuit Diagram of PWM based Matrix 
Converter 


Figure 6: Matrix Converter 



Figure 7: SVM based Matrix Converter 


4.1 Case Studies 

4.1.1 Test Case I A: Fixed Input Frequency 
and Variable Output Frequency 


In the test case I A, input frequency is fixed and out- 
put frequency is varied. It is observed that change in 
output frequency does not affect the input frequency 
when PWM and SVM based AC to AC matrix con- 
verters are employed as shown in figures 10, 11 and 


12 . 


4.1.2 Test Case IB: Variable Input Frequency 
and Fixed Output Frequency 

In the test case IB, output frequency is fixed and input 
frequency is varied. It is observed that change in input 
frequency does not affect the output frequency when 
PWM and SVM based AC to AC matrix converters 
are employed as shown in figures 13, 14 and 15. 

4.1.3 Test Case IIA: Fixed Input Voltage and 
Variable Output Voltage 

In the test case III A, input voltage is fixed and out- 
put voltage is varied. It is observed that change in 
output voltage does not affect the input voltage when 
PWM and SVM based AC to AC matrix converters 
are employed as shown in figures 16, 17 and 18. 

4.1.4 Test Case IIIB: Fixed Output Voltage 
and Variable Input Voltage 

In the test case III, output voltage is fixed and in- 
put voltage is varied. It is observed that change in 
input voltage does not affect the output voltage when 
PWM and SVM based AC to AC matrix converters 
are employed as shown in figures 19, 20 and 21. 
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Figure 9: Circuit Diagram of SVM based Matrix Con- 
verter 



Figure 10: 20Hz Output Signal with 50Hz Input Sig- 
nal Frequency. 



Figure 11: 80Hz Output Signal with 50Hz Input Sig- 
nal Frequency. 



Figure 12: 100Hz Output Signal with 50Hz Input Sig- 
nal Frequency. 



Figure 13: 20Hz Input Signal with 50Hz Output Sig- 
nal Frequency. 



Figure 14: 80Hz Input Signal with 50Hz Output Sig- 
nal Frequency. 


4.1.5 Test Case III A: FFT Analysis for Pulse 
Width Modulation Based Matrix Con- 
verter with Fixed Input Frequency 


Fast Fourier Transform Analysis for PWM based ma- 
trix converter is carried out as shown in figures 22, 
23 and 24. It is found that this matrix converter 
maximizes the fundamental content of the character- 
istics by reducing the harmonic content. Any desired 
output frequency with maximum fundamental content 
can be obtained with this model by fixing the supply 
frequency. 


4.1.6 Test Case III B: FFT Analysis for Pulse 
Width Modulation Based Matrix Con- 
verter with Fixed Output Frequency 


Fast Fourier Transform Analysis for PWM based ma- 
trix converter is carried out as shown in figures 25, 
26 and 27. It is found that this matrix converter 
maximizes the fundamental content of the character- 
istics by reducing the harmonic content. Any desired 
input frequency with maximum fundamental content 
can be supplied with this model by fixing the output 


51 


Automatic Control and System Engineering Journal, ISSN 1687-4811, Volume 16, Issue 2, ICGST, Delaware, USA, December 2016 



Figure 15: 100Hz Input Signal with 50Hz Output Sig- 
nal Frequency. 



Figure 19: 110V Input Signal with 220V Output Sig- 
nal Frequency. 



Figure 16: 110V Output Signal with 220V Input Sig- 
nal Frequency. 



Figure IT: 440V Output Signal with 220V Input Sig- 
nal Frequency. 



Figure 18: 1000 V Output Signal with 220V Input Sig- 
nal Frequency. 



Figure 21: 1000V Output Signal with 220V Input Sig- 
nal Frequency. 



Figure 22: FFT analysis of 20Hz Output Signal with 
50Hz Input Signal. 


frequency or in other words, fixed frequency applica- 
tions can work smoothly irrespective of the variation 
in the supply frequency by deploying this model. 
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Figure 23: FFT analysis of 80Hz Output Signal with Figure 25: FFT analysis of Fixed Output 50Hz with 
50Hz Input Signal. 20Hz Input Signal. 
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Figure 24: FFT analysis of 100Hz Output Signal with 
50Hz Input Signal. 
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Figure 26: FFT analysis of Fixed Output 50Hz with 
80Hz input Signal. 


4.1.7 Test Case IV A: FFT Analysis for Space 
Vector Modulation Based Matrix Con- 
verter with Fixed Input Frequency 


Fast Fourier Transform Analysis for SVM based ma- 
trix converter is carried out as shown in figures 28, 
29 and 30. It is found that this matrix converter 
maximizes the fundamental content of the character- 
istics by reducing the harmonic content. Any desired 
output frequency with maximum fundamental content 
can be supplied with this model by fixing the input 
frequency. Total harmonic distortion is lower as com- 
pared to PWM based matrix converter. 


4.1.8 Test Case IV B: FFT Analysis for Space 
Vector Modulation with Fixed Output 
Frequency 

Fast Fourier Transform Analysis for PWM based ma- 
trix converter is carried out as shown in figures 31, 
32 and 33. It is found that this matrix converter 
maximizes the fundamental content of the character- 
istics by reducing the harmonic content. Any desired 
input frequency with maximum fundamental content 
can be supplied with this model by fixing the output 
frequency or in other words, fixed frequency applica- 
tions can work smoothly irrespective of the variation 
in the supply frequency by deploying this model. To- 
tal harmonic distortion is lower as compared to PWM 
based matrix converter. 
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Figure 27: FFT analysis of Fixed Output 50Hz with 
100Hz input Signal. 



Figure 28: FFT analysis of 20Hz Output with 50Hz 
Input Signal. 



Figure 29: FFT analysis of 80Hz Output with 50Hz 
Input Signal. 


5 Conclusion 


In this research, PWM and SVM based three phase 
AC to Ac matrix converters are developed. Simula- 
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Figure 30: FFT analysis of 100Hz Output with 50Hz 
Input Signal. 
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Figure 31: FFT analysis of output with 50Hz Output 
and 20Hz Input Signal. 


tions are performed on Matlab. It is found out that 
PWM and SVM based matrix converter can be de- 
ployed to achieve any desired output and input char- 
acteristics. These converter are highly applicable for 
adjustable speed drives or variable frequency drives 
because in these proposed models it is shown that ir- 
respective of any supply frequency, variable desired 
frequency can be attained. Secondly, this proposed 
model has revolutionized the applications requiring 
variable voltages. Any desired output voltage can be 
achieved without taking into consideration of input 
voltage and vice versa. Interconnected systems can 
be free of synchronization issues by employing these 
converters. Further, harmonic contents are greatly re- 
duced and it has maximized the fundamental content 
of the desired characteristics. SVM based matrix con- 
verter offer low total harmonic distortion than PWM 
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Figure 32: FFT analysis of output with 50Hz Output 
and 80Hz Input Signal. 




Splnrfflri. 1 rjfcln*. rFT winrlpw (-n rwfy I rfrUt 


-JO 

-K 

vA/' i 



6.01 603 00J OCW 0.04 6 

Hrti*W 

06 


Tmtnm* £QMl) . m 1 MiE* 


16 

Ifl ■ 

TV • 

1 0 ■ 

5 * 

L 



„ , 

in am an m sn eo m iwo an ip 

r int|MBncT CHjJ 

an 


Figure 33: FFT analysis of output with 50Hz Output 
with 100Hz Input Signal. 
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Abstract 

Switched reluctance motor (SRM) has outstanding 
advantages such as simple structure, absence of 
commutator, high torque. However, torque ripple and 
loud noise are the biggest disadvantages of this SRM. 
Thus, applications in some special areas in industry are 
restricted. To reduce torque ripple and improve output 
power, this article provides solutions for control and 
change the structure of SRM. At the same time, it also 
proposes a new solution: using of new materials - 
amorphous materials in manufacturing rotor and stator 
instead of normal Silicon steel materials. With special 
physical properties, the amorphous materials used in 
motor manufacture give a greater speed and lower loss. 

Keywords switched reluctance motor, tor ue ripple, 
amorphous materials. 

Nomenclature 

AMN Associative Memory Networks 

APC Angular Position Control 

CCC Current Chopping Control 

CVC Voltage Chopping Control 

FEM Finite Element Method 

PI Proportional Integral 

PID Proportional Integral Derivative 

PWM Pulse Width Modulation 

RBF Resting Bitch Face 

SMC Sliding Mode Control 

SRM Switched Reluctance Motor 

TSF Torque Share Function 

1. Introduction 

Switched reluctance motor (SRM), a type of simple 
structure motor, was proposed since 1842. However, 
SRM operates under the principle of switching electric 
current for each phase, so it was very difficult to 
control SRM at that time. By the 80s, with the rapid 
development of power electronics, microprocessors and 
control technology, SRM started becoming the choice 
for high speed control systems. 



Switched reluctance motor has the advantages such as: 
simple structure, high reliability, the motor can work at 
an extremely fast speed. Thanks to permanent magnet 
and winding are absent in the rotor, the rotor 
temperature can be higher than other types of motor. 
Torque direction doesn’t depend on current direction, 
so can simplify the inverter, reduce the cost of the 
system. As there isn’t overcurrent that can damage the 
capacity valve in SRM, the inverter has a high 
reliability. Big starting torque, high speed controlling 
performance, there is no impact of the current in the 
induction motor at the time of starting. Magnetic circuit 
of motor can work in the linear region and the 
saturation of the B-H magnetic characteristic curve, 
using the maximum capabilities ferromagnetic 
materials, power ratio with weight is large. With the 
above advantages, SRM is widely used in industry and 
some specific sectors such as manufactured generators 
for aircraft launch, the centrifuges require high speed, 
dynamic the submarine base with large startup torque, 
etc. 

However, due to salient pole structure of rotor such 
motor flux has strong non-linearity, torque ripple is 
large and highly noise. Therefore, in some cases, the 
applications of SRM are limited. In order to overcome 
these restrictions of torque ripple and improve the 
quality of motor, based on analysis, comparison and 
summarization of domestic and international relating 
documents, the authors of this article have figured out 
some solutions to reduce torque ripple including 
solutions of controlling and solutions of changing 
motor structure. At the same time, the authors also 
proposed a new solution - using new materials in the 
motor manufacturing to improve the characteristics of 
this switched reluctance motor. 

2. Mathematical model of switched 
reluctance motor 
2.1. Voltage equation 

Survey of SRM operating status requires a 
mathematical model. Since SRM is a non-linear 
system, it demands a suitable model which represents 
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the dynamic process in various operating conditions. 
Assuming the phase of the motor including: resistance 
of winding R, inductance L(i, 0), the voltage u supply 
to the phases by total voltage drop and the rate of flux 
\|/(i, 0), effects of mutual inductances can be neglected, 
Faraday's law gives: 


dy(i,0) _ i 

dt 


<r+U =Ri + 


dy/(\,6) 

dt 


where = (/,#)/ 


( 1 ) 

( 2 ) 


Switched reluctance motor structure is symmetrical so 
flux circulation cycle 27i/nR (nR is the number of rotor 
poles). Hence, there are two concepts to be 
distinguished: mechanical angle 0 and electrical angle 
(n R 0). Mechanical angle indicates the rotation angle of 
the rotor position, electrical angle shows angle 
calculating phase current and voltage. Voltage equation 
can be rewritten in the form of electrical angle as 
follows: 


U = Ri + 


dyjnAi) 

dt 


(3) 


Consider the three phase motor, each phase mismatches 

2 jz 

of a electrical angle: 6. = . Continue to ignore 

3 n R 

mutual inductance, equations for the phase voltages are 
as follows: 


u a = Ri *+-^v'( n R0,i s ) 
<u b =Ri b +j tV (n K (0-0 s ),i b ) (4) 

u c=Ric+j t V / ( n « (0-20M 


2.2. Torque equation 

In any rotor position of a phase, according to the 
principle of total of magnetic energy, the torque 
formula of the SRM is calculated as follows: 


do 

0 is rotor azimuth 

1 is phase current 

W (0, i) is the total of magnetic energy 


(5) 


Total of magnetic energy is defined as the area under 
the curve of magnetization shown in Figure 1, it can be 
expressed as follows: 

i 

W'(0,i) = jy/(0,i)di (6) 

0 



Figure 1. The chart outlines the total magnetic 


In certain cases, SRM is not affected by magnetic 
saturation and inductance is independent of electric 
current, depends only on rotor position (2), the torque 
generated as follows: 

1 i 2 

T(0,i) = — — f (0,i)di 2 (7) 

2d0{ 

If L changes linearly with rotor position, torque 
equation is as follows: 

rAr-LM (8) 

2 dt 


Torque of three phase motor equal to total of phases 
torque: T e =T a +T b +T C . 

2.3. Mechanical equation 

According to the torque balance theorem of motor: 

T e~l=J % (9) 


T e is total torque 
Ti is load torque 
flm is speed 
J is inertia torque 
B is constant friction 

Based on three equations above, we can build model of 
switched reluctance motor. This model based on the 
change of voltage equations: 


%fV(n R 0,i a ) 

5 'a 

di ~ R 'o " JZ ~ '» a + ' u » 

u b C0 

dt ^r(n„(0-0 s ),i b ) 

dl b 

di c ~ R ’ C -^ 0 V / ( n K (0- 2 0 s )Jc) n RV+ u c 

dt ^- v (n R (0-20 s ),i c ) 

dl c 


( 10 ) 
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3. Solutions to improve the characteristics 
of switched reluctance motor 

Switched reluctance motor torque ripple is really 
complicated and affected by many factors. Therefore, 
to reduce torque ripple, strategic solutions to be given 
are using SRM control method with definite motor 
structure and parameters [1] or changing motor 
structure. In this paper, the authors have presented, 
analyzed and compared with domestic and international 
researches relating to reduction of torque ripple and 
improvement of SRM characteristics. 

3.1. Solution using control method of switched 
reluctance motor 

According to equation (1) and (8), motor control 
method can be performed by current control i, or 
adjust angle 0. These are two important parameters to 
calculate the voltage and torque of SRM. 

3.1.1. Traditional control strategy 

There are three ways in typical traditional control 
including angular position control (APC), current 
chopping control (CCC) and voltage chopping control 
(CVC). 

The angular position of the control means at a certain 
voltage across the winding provided by changing 0 on 
and 0 O ff of main switch to change power-on and power- 
off. Thus, to control closed loop speed, it needs to 
adjust the phase current waveform signal. Opening 
angle 0 on and off angle 0 O ff can be revised, however in 
fact, we often have to fix 0 O ff and change 0 on according 
to control mode. 0 on and 0 O ff must be adjusted 
separately, then SRM is able to achieve optimum 
operating angle position control (APC) mode. There 
are some advantages of this control method such as 
greater torque adjustment range, multi phase windings 
simultaneously powered, high motor efficiency, 
suitable for higher speeds but unsuitable for low speed 
operation. 

Current chopping control is mainly used for motor 
running at low speed. Because of increasing rapidly of 
the current, in order to avoid damage caused by current 
pulse, it requires to limit the peak current. This control 
method rarely involves 0 on and 0 O ff control, so select the 
position of each phase winding's conduction to control 
the current by chopper. The main advantages of this 
control method are simplicity and accuracy, reduction 
of torque ripple, suitable for the low speed operation of 
fixed load motor. 

Voltage chopping control means keeping 0 on unchanged 
and changing 0 O ff, the power switching devices operate 
in pulse width modulation (PWM) mode. Keep the 
pulse period T fixed, simply adjust the duty cycle of the 
PWM waveform, thereby adjust the size of the voltage 
across the winding, which cause changes in winding 
phase current to achieve the regulation of motor speed. 
Advantages of the method are followings suitable for 
both high speed and low speed operation, enable quick 
adaptation to change of the load. 

3.1.2. Torque distribution strategy 

Torque distribution strategy shows that by defining the 
torque share function (TSF), electromagnetic torque 
corresponding to the current phase can be distributed 


reasonably to ensure each phase torque being at a 
constant value. Accordingly, phase current values can 
be calculated and combined with appropriate control 
strategies to reduce torque ripple. 

Reference [2] proposes an optimized solution based on 
conventional torque distribution strategy to maximize 
the torque/ampere ratio within the allowable error 
range of torque ripple in adjusting the opening angle 0 on 
and off angle 0 O ff of SRM. Under the premise of 
guaranteed output torque, to achieve the optimal 
purpose, it is important to ensure lower the winding 
phase current. After testing, the results show that the 
optimized ratio of torque and current will improve the 
system performance. 

For the SRM torque ripple, Reference [3] adopts the 
torque distribution based on controller to obtain the 
desired torque of each phase, and then synthesize the 
constant torque through each phase, effectively 
eliminate torque ripple. 

Reference [4] combines the current torque distribution 
and current to achieve efficient suppression of the 
torque ripple. Up to a point, the method has achieved 
certain results, but choosing to divide torque has 
omitted necessary physical factors. 

In order to extremely minimize the torque faults, 
reference [5] dissects the optimization by using flux 
linkage as optimized variables. With discrete method, 
the torque reducing problem is transferred into a 
mathematical programming issue. Its feasibility has 
been proved by simulation. But this method ignores 
the correlation between electromagnetic and winding. 

3.1.3. Intelligent control 

Intelligent Control is a type of control of which object 
is mathematical non-linear, its properties are strong 
self-learning ability and adaptive capacity. Applying 
the intelligent control theory in designing SRM 
controller is one of the effective ways to solve thorny 
issue of SRM meaning torque ripple reduction. 
Performance control of SRM is significantly improved 
a. Fuzzy control 

Typical application of fuzzy control is an adaptive 
fuzzy control strategy of minimum SRM torque ripple 
proposed by Sayeed [6]. 



Figure 2. The block diagram of adaptive fuzzy torque 
control 

Figure 2 shows the adaptive fuzzy SRM torque control 
system regarding position as an input, phase current as 
an output. Changing the membership function gave 
broken phases in the most appropriate area. The initial 
value of the fuzzy parameter is randomly selected. 
During the operation, the parameters are adjusted by 
the real time controller, able to adapt to the motor 
characteristics of any changes. 
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Reference [7] applies fuzzy neural algorithm constantly 
adjusting the size of the current compensation signal to 
optimize the phase current waveform and then to 
suppress torque ripple. Because of the difficulty in 
measuring the real time of the dynamic torque and 
expensive torque sensor, it is hard to apply this method. 
Reference [8] through the fuzzy neural network speed 
controller determines the desired torque output, 
calculates the current value of each phase. 

Reference [9] proposes a new method of mathematical 
modeling of a motor switched reluctance, based on the 
static inductance curve and torque - angle characteristic 
curve, apply the adaptive network based on fuzzy 
inference system to the whole modeling of SRM. 
b. Neural network 

Reducing torque ripple based on neural network was 
first proposed by JGO Donovan [10]. SRM non-linear 
characteristics, SRM torque - current - angle 
relationship are analyzed based on the knowledges, 
experiences of professional experts. 

Reference [11] applies AMN (associative memory 
networks) neural network online to get the motor phase 
current waveform optimization. Disadvantages of this 
approach are expensive torque sensor, inaccurate 
measurements at high speeds and no torque feedback 
during normal operation. 

Reference [12] proposes combination of traditional PID 
with neural network to make adaptive PID control 
applied for serious non-linear SRM. The non-linear 
transformation properties and parallel processing 
ability at a high degree of neural network made it 
suitable for the establishment of non-linear predictive 
model parameters to predict. Through predicting the 
parameters of the control system, the dynamic response 
performance of the system can be improved as shown 
in Figure 3. The results show that the system is stable, 
quick, accurate, strong anti-disturbance and easy 
control. 



Figure 3. Block diagram of the non-linear model based on 
RBF neural network 

Reference [13] proposes a SRM instantaneous torque 
control method which is based on radial basis function 
neural network. Using data generated from the SRM 
dynamic simulation for RBF neural network's offline 
training, learning to optimize current waveforms under 
different speed and torque, and then RBF network is 
used for controlling motor torque. Experimental results 
show that the control strategy can effectively reduce 
SRM torque ripple with high control accuracy, able to 
adapt to speed changes. 


3.1.4. Sliding mode control 

Reference [14] proposes a new method to minimize 
torque ripple of a switched reluctance motor, which is 
based on the control of the sum of the square of the 
phase currents by using only two current sensor and 
analog multipliers. In addition, the sliding mode control 
(SMC) technique has been applied to SRM speed 
control loop that compensates the low frequency 
oscillations on the torque output. The results show that 
the continuous sliding mode controller is effective 
compared to PI or Fuzzy controllers in reducing the 
torque ripple of the motor, compensating for non-linear 
torque characteristics and making drive insensitive to 
parameter variations as well, low frequency vibration in 
torque variation has been eliminated. 



Figure 4. Block diagram for the speed control with SMC 


In addition to the methods described above, the issue of 
SRM control to improve the work quality, reduce 
torque ripple which researchers overcome by the other 
control methods such as linearization control, Iterative 
Learning Control, Internal Model Control, etc. 

3.2. Solution for change of switched reluctance 
motor structure 

Equation (5) and (9) show the torque of motor directly 
affected by the magnetic flux density which is 

closely related to the motor structure. Therefore, motor 
structural change solution is presented to improve 
torque characteristics and SRM performance. 

3.2.1. Change the shape and size of stator and rotor 
pole 

General purposes of designing switched reluctance 
motor is to maximize torque, reduce torque ripple by 
optimizing control manner and strategy. In terms of 
trade, however, the design target is not only to meet the 
requirements of torque at both low and high speed but 
also to reduce the cost of motor production. Reference 
[15] uses the finite element analysis method (FEM) to 
simulate the shape of SRM with many different 
configurations by changing the shape of stator and rotor 
poles. 

Rotor is usually made of ferromagnetic materials 
(silicon). Ferromagnetic materials have a direct 
influence on making flux. The author changed iron 
ratio in the rotor (62%, 76%, 91%). In there, this rate of 
SRM rotor is 62% for max flux (0.03 lWb). Thus, 
because of strong nonlinearities, inductance dependent 
of angle 0 and current i, changing ferromagnetic 
materials rate in SRM rotor contributes only to 
optimize calculations in design while has not shown a 
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clear positive correlation between ferromagnetic 
materials rate and flux in the motor. 

The resizing rotor and stator poles size related directly 
to the change in order to achieve the ratio y (y is the 
ratio between rotor pole angle and pole pitch) and P (P 
is the ratio between stator pole angle and pole pitch). 
The stator and rotor pole angles are important variables 
in SRM design, it affects to the motor torque and it also 
has been studied to find precisely range of the stator 
and rotor pole angle in practice. A parameter 
obtainedfrom generation ofresults in changing the poles 
is the ratio between the polar angle and pole pitch (pole 
pitch is the angle between two consecutive stator poles 
or rotors). Survey results in 3 phase SRM, reference 
[15], give the optimal value y and P (y = 0:38, P = 0.5) 
for the highest performance. 

The advantage of this method is figured out the optimal 
parameter y, P and torque increased. However, the 
torque also depends on many other motor configuration 
parameters; motor is only surveyed to stop at 3 phase 
SRM 6/4 poles and 8/6 poles, low speed range. 

3.2.2. Change rotor and stator pole surface 

The non-linear torque in the overlap region is the main 
cause of torque ripple. Also, the radial force is 
generated in the whole overlap area of stator and rotor 
poles and the tangential force is generated at the edges 
of stator and rotor poles, respectively. To reduce the 
fringing flux at the overlap angles, a rotor shape made 
by method of laminated layer of the asymmetry rotor is 
proposed [16], as Figure 4a. To have the asymmetry 
angles, the width of each steel sheet at rotor pole arc 
needs to be different. Research results conducted with 
SRM 12/8 poles show that the torque ripple was 
reduced by 15% against the initial SRM and increased 
by 12.5% compared with normal structures. However, 
switched reluctance motor was tested with quite small 
capacity, just 1.5kW. 




Figure 4. Surface structure of rotor pole (a) and stator 
pole (b) 


Stator pole shape has been changed in the reference 
[17]. Stator pole surface is constructed in tip shape as 
shown in Figure 4b, flared wide at the tip to increase 
the amount of magnetic flux passing, the magnetic flux 
density exists in the stator poles increased, resulting in 
increased torque. The author used multi domain 
simulation environment (ANSYS workbench) to 
simulate the magnetic field, torque of switched 
reluctance motor. 

3.2.3. Change the number of rotor pole 

The number of stator pole in initial SRM is more than 
number of rotor pole (6/4, 8/6, 12/8, etc). However, 
6/10 SRM introduced in reference [18] is a new motor 
concept, in which the rotor poles are more than the 
stator poles, focusing on creation of greater torque, 
improving torque density, reducing torque ripple while 
maintaining continuous thermal energy density. These 
advantages are thanks to structure of 6/10 SRM with 
pole pitch and rotor angle smaller than traditional 
structure. The more aligned positions between rotor and 
stator are created the more flux with maximum value 
increased. 

The author has made a comparison of 6/10 SRM and 
6/4 SRM about the flux in stator winding, thermal 
efficiency, torque density and torque ripple. With 6/10 
SRM, torque performance increases (up to 119.82% 
compared with the 6/4 SRM). However, the winding 
temperature is more than 40% and torque ripple is not 
reduced. 

Reference [19] showes the new structure of SRM with 
equal number of rotor poles and stator poles (6/6 pole). 
However, with such a structure of pole, the motor can 
not be started by itself. But this is also a special point 
because this motor uses only single phase. Thus, the 
system has a simple converter having two important 
roles: to improve the input power factor and to 
stimulate the motor phase. A test was conducted for 6/6 
pole motor with rated power of 400 W and speed of 
1800 rpm. However, this study is just applied for the 
survey of current and voltage during engine operation 
with the expectation of high performance at low cost 
of system. 

3.2.4. Change the number of stator teeth 

In principle, it is expected that increasing the number of 
teeth per stator pole may be increase more advantages, 
limit disadvantages and enlarge widely certain 
applications of switched reluctance motor. It seems 
generated torque has increased owing to increasing 
number of teeth, the losses in the motor are limited, 
particularly at high speeds. These losses are lower than 
that in SRM with only one tooth per stator pole. The 
first SRM with many teeth per stator pole which is 6 
pole 3 phase 12/10 slotted motor with 2 teeth per pole 
stator. Other version of 4 teeth per stator pole is 3 phase 
motor with 22 tooth rotor and 24 stator teeth [20], as 
shown in Figure 5. Both these motors are designed with 
same frame size and air gap for easy comparison. The 
simulation result in two bands speed is 300 rpm and 
750 rpm showing greater motor performance and high 
torque. However, at a speed of 300 rpm, 24/22 SRM 
has higher performance than 12/10 SRM; at speed of 
750 rpm, the result is opposite. 
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Besides the methods described above, some researches 
continue to be made to improve the characteristics of 
SRM as adjusted inclined slotted, perforated on the 
rotor teeth, tilted stator pole, rotor, etc. Those methods 
have contributed to improve the work quality of the 
switched reluctance motor in certain cases, but the 
contribution in reducing torque ripple is not thoroughly 
made. 



Figure 5. Stator pole structure many teeth 


3.3. Propose a solution using new materials 

Equation (2) shows the relationship between magnetic 
flux \| / with current i and angle 0. This flux curve is 
non-linear, depends on the materials of rotor, stator. In 
equation (5), (6), the effect of magnetization of 
materials to torque, torque ripple is quite clear. 
Therefore, changing the motor manufacturing materials 
is an effective solution in order to improve motor 
characteristics. 



Figure 6. Rotor structure used aluminum materials and 
magnetic field lines 


Currently, both switched reluctance motor in particular 
and other motor types in general use electrical steel in 
manufacturing rotor and stator. Electrical steel is a soft 
ferromagnetic, large magnetic permeability, small 
degaussing magnetic, low hysteresis loss (due to 
narrow magnetic hysteresis curve). However, 
ferromagnetic materials exists limitations such as low 
operating frequencies (silicon steel cores can be used 
only for low frequency, at high frequency, losses will 
be greater due to low resistivity of silicon steel), small 
coercivity force, etc. Moreover, due to the requirements 


of environmental protection and saving energy, 
reducing core losses and magnetic losses, the solutions 
using other materials instead of ferromagnetic materials 
started to be proposed. 

Reference [21] gives a new type of SRM model used 
unibody aluminum for rotor core, the purpose is to 
reduce motor weight, as shown in Figure 6. At the same 
time, to create magnetism in the rotor, the magnetic 
core is embedded in the rotor blocks. This design has 
been used to make the low speed motor. The simulation 
results shows that torque did not significantly increase, 
due to the fact that the magnetic field lines didn’t go 
through the central rotor (rotor core made of aluminum) 
causing the low density of magnetic flux. However, 
thanks to removing the radial magnetic force, noise is 
significantly reduced. 

Based on the results of research using new materials 
being low-iron-loss such as Super E-Core brand, 
reference [22] also made comment that strip steel has 
improved performance of the SRM and prevent steel 
loss. 




Applied Field | A m] 

Firgue 8. Saturation of flux density with the different 
materials 


Through the synthesis of researches on some new 
materials with good magnetic and high conductivity, 
authors have proposed solutions using new materials in 
order to replace conventional ferromagnetic material. 
This new material is amorphous alloys. Amorphous 
alloy is a kind of alloy made of iron or cobalt (Co), in 
the amorphous state, so its resistivity is much higher 
than crystalline alloy. At the same time, this material is 
resistant to corrosion, high mechanical strength, can be 
used at higher operating frequencies than crystal 
material of metal base. Without crystal structure of 
amorphous material permits elimination of, anisotropic 
so its properties is highly soft magnetic. Cobalt 
amorphous materials has magnetostrictive by 0 so it has 
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a very small coercive force, can be compared with 
nanocrystalline materials. Magnetic hysteresis curve of 
this material is much narrower than silicon steel - soft 
magnetic materials (Figure 7). Natural thickness of the 
amorphous steel sheets is very small around ~ 0.03 mm 
compared with ~ 0.3 - 0.5 mm of silicon steel, 
resistivity greater - 130- 170 juQcm compared with ~ 
50 - 60 /jQcm of silicon steel. Therefore, steel core 
loss reduced significantly compared with silicon steel. 
To demonstrate the effectiveness of the research 
results, the research team has developed the 
characteristic curves for evaluation that is described in 
figure 8. Thus, it is shown that amorphous steel is 
capable of spontaneous magnetization and long-lasting 
magnetizing energy which increased the magnetic flux 
density and improved motor torque (Figure 8). 
Moreover, amorphous material weight is very light 
useful for motor operation at much higher speed than 
normal motor. Analytical result shows that steel core 
losses in amorphous materials decreased by about 20% 
compared with silicon steel, magnetic flux density 
increased by over 12% and the motor speed can be up 
to 140% compared to the motor manufactured by 
conventional silicon steel. The analytical result is 
shown in Table 1. 


Table 1. Result of the study compared with the 
traditional SRM [23] at speed n = 3000 rpm and current 
1 = 4 A. 


Motor specs 

New 

Traditional 

SRM 

SRM 

Copper losses (W) 

2.24 

2.8 

Core losses (W) 

1.2 

1.4 

Motor efficiency (%) 

83.7 

75.4 

Aligned inductance (mH) 

4.2 

3.75 

Unaligned inductance (mH) 

0.78 

0.7 

Average torque (Nm) 

0.053 

0.047 


Therefore, the research and design of switched 
reluctance motor using amorphous materials is an 
important and necessary solution. However, it is 
necessary to evaluate the relationship of the magnetic 
field in the motor to give a specific design technology 
and manufacture using amorphous materials. Research 
result announced by the authors in this paper is initial, 
the standardized assessment results through the 
simulation of physical transformation process of the 
motor parameters will be announced in the appropriate 
time. 

4. Conclusion 

This article presents different solutions in order to 
improve the characteristics of switched reluctance 
motor. With the development of various intelligent 
control theories, SRM controller has achieved good 
results in reducing torque ripple. At the same time, 
research of solution changing the internal structure of 
motor which will be continued because of its effective 


solution to strictly overcome the disadvantages of the 
switched reluctance motor. 
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